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Theoretical description of reactions in general, and the theory for (d,p) reactions, in particular, 
needs to advance into the new century. Here deuteron stripping processes off a target nucleus 
consisting of A nucleons are treated within the framework of the few-body integral equations theory. 
The generalized Faddeev equations in the AGS form, which take into account the target excitations, 
with realistic optical potentials provide the most advanced and complete description of the deuteron 
stripping. The main problem in practical application of such equations is the screening of the 
Coulomb potential, which works only for light nuclei. In this paper we present a new formulation 
of the Faddeev equations in the AGS form taking into account the target excitations with explicit 
inclusion of the Coulomb interaction. By projecting the (A + 2)-body operators onto target states, 
matrix three-body integral equations are derived which allow for the incorporation of the excited 
states of the target nucleons. Using the explicit equations for the partial Coulomb scattering wave 
functions in the momentum space we present the AGS equations in the Coulomb distorted wave 
representation without screening procedure. We also use the explicit expression for the off-shell 
two-body Coulomb scattering T-matrix which is needed to calculate the effective potentials in the 
AGS equations. The integrals containing the off-shell Coulomb T-matrix are regularized to make 
the obtained equations suitable for calculations. For AW and nucleon-target nuclear interactions 
we assume the separable potentials what significantly simplifies solution of the AGS equations. 

PACS numbers: 21.45.-v,24.10 Eq, 24.50 +g, 25.45 -z 
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I. INTRODUCTION 

Scattering of composite projectiles such as the deuteron off a composite target nucleus can in principle be described 
by the N-body scattering theories (see, e.g., Ref. [l|), which allow the simultaneous treatment of all reaction channels 
(including partial and complete breakup) in an exact, unique, N-particle-unitarity-preserving manner. These N-body 
equations can be reduced to two-cluster equations, as proposed for the first time in Ref. Q, the only input necessary 
being the elementary two-body transition amplitudes. However, due to the inherent complexity of the resulting 
equations it is always desirable to take resort to some kind of simplifying approach. For reactions in which, besides 
two-body, also three-body channels are known to be important a promising candidate is the three-body Faddeev 
integral equations theory written in the Alt-Grassberger-Sandhas (AGS) form 0, Q. But even the three-body AGS 
equations are so complicated that a more simplified CDCC approach became quite popular in the last 30 years (see 
[3-0] and references therein). Of course, this approach does no longer constitute a rigorous theory, the approximative 
character becoming the more apparent the more important the internal degrees of freedom are which can be excited 
in the energy domain under consideration. 

If the internal structure of the target can no longer be neglected, in order to still be able to work with the manageable 
three-body theory, the possibility of excitation must - at least approximately - be taken into account. For the first time 
it has been done in Q . Deuteron stripping processes off a target nucleus consisting of A nucleons were treated within 
the framework of the generalized few-body Faddeev integral equations theory written in the AGS form. Generalization 
of the AGS equations is achieved by taking into account the excitation of the target. To reduce the (A + 2)-particle 
problem to the much simpler three-body problem all the operators acting in the (A + 2)-particle space were projected 
onto the three-particle space. Obtained generalized Faddeev equations couple all rearrangement (reaction), inelastic 
and elastic amplitudes. The transition amplitudes for all interesting three-body processes were obtained, whether the 
target nucleus is in its ground or in some excited state before and/or after the collision. The practical application 
was done for the deuteron stripping reaction 12 C(d,p) 13 C at the deuteron bombarding energies of 4.66, 15 and 56 
MeV. The two-body-type T-operators for the nucleon-nucleus subsystem were calculated from multichannel equations 
to account for the excitation and de-excitation of the nucleus in nucleon-nucleus scattering. When inserted in the 
three-body-type integral equations this feature is automatically introduced also into the three-body dynamics. The 
complexity of the AGS equations is significantly reduced by using separable potential approach for NN and NA 
interactions. However, the Coulomb p — A interaction was neglected when solving the AGS equations and was taken 
into account only approximately by multiplying each calculated purely nuclear, partial wave reaction amplitudes by 
the initial- and final-state Coulomb distortion factors before summing up the partial wave series. Definitely, in such an 
approach we neglect the Coulomb-modified vertex form factors describing the subsystem p— A and effective potentials 
in the AGS approach described by the triangular diagrams containing the Coulomb p — A scattering amplitude as the 
four-ray vertex. 

After our work Q a significant advance in the application of the Faddeev three-body approach was achieved in 
works l9 Hllll. where AGS equations for stripping reactions on different targets were solved with realistic potentials. 
In works [Ha. [l3j the AGS calculations were compared with the adiabatic distorted wave approach (ADWA) for 
deuteron stripping reactions to estimate the accuracy of the conventional ADWA, which provides a simplified but 
practical version of the CDCC [l4| • The advantage of these works compared to our work Q was usage of the realistic 
potentials. Also the Coulomb interaction was included using the screening procedure, which has been applied earlier 
in our work for p+d scattering [l5j . Inclusion of the screening procedure requires higher screening radii when charge of 
the target increases to get the convergence. That is why the application of the Coulomb potential screening procedure 
was successful only for targets with charge Z < 20. Besides, the screening procedure cannot be always a reliable 
remedy to solve problems with charged particles, because limiting the range of the potential may lead to the loss of 
information about the very nature of the field creating the Coulomb potential. Neglect of the internal structure of 
the target is another setback in the AGS calculations in f9l-[l3l]. 

In this paper we present a new formulation of the generalized Faddeev equations in the AGS form for the deuteron 
stripping, which includes explicitly the Coulomb interactions and target excitations. The Coulomb interaction in 
the AGS approach appears in the three-ray vertex form factors in the effective potentials and in the four-ray vertex 
in the triangular diagrams. Applying the two-potential equation allows us to remove the non-compact singularity 
in the triangular diagram describing the elastic scattering and containing the p — A Coulomb scattering amplitude. 
Besides the AGS equations can be rewritten in the Coulomb distorted wave representation (l6l - [l8j . in which the 
reaction amplitudes and the effective potentials are sandwiched by the Coulomb distorted waves in the initial and 
final states. Applying the regularization procedure we obtain the expression for the effective potentials in the AGS 
equations, which are free of the singularities caused by the Coulomb distortions in the initial and final states. We 
also investigate the off-shell Coulomb scattering amplitude and show that the Coulomb-modified form factors in the 
transfer amplitudes and in the triangular diagrams don't contain non-integrable singularities. The target excitation is 
taken into account following the formalism developed in The final generalized matrix AGS equations are written 
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in the form which includes explicitly Coulomb interactions, target excitations and spins of the particles. Because 
the solution of the AGS equations is greatly simplified for separable potentials, we use separable potential approach 
assuming that the adopted separable potentials will approximate realistic potentials NN and NA potentials. The 
calculations of the obtained equations and comparison with the experimental data will be presented in the following 
up papers. 

It is important to underscore that AGS equations in the Coulomb distorted wave representation have been derived 
for the first time. Previously possibility of derivation of these equations in the Coulomb distorted wave representation 
was considered in fl6l [l8j , but obtained equations were not genuine integral equations because the amplitudes under 
the integral and on the left-hand side were different functions. Here, using a different off-shell continuation we derive 
genuine AGS equations in the Coulomb distorted wave representation. Another important topic of our research is 
investigation of the singularities of the integrals containing the off-shell Coulomb scattering amplitude. The com- 
pactness of the AGS equations for charged particles with repulsive interactions was proved in However, the 
practical application of this result requires regularization of the integrals containing the off-shell Coulomb scatter- 
ing amplitudes. Note that in the case under consideration, when only two particles are charged, only one off-shell 
Coulomb scattering amplitude of the proton-nucleus scattering is needed. Regularization of the integrals containing 
the off-shell Coulomb scattering amplitude is also important because in a standard procedure involving the screen- 
ing of the Coulomb potential is tacitly assumed that, in the limit of the screening radius R — > oo, all the integrals 
containing Coulomb scattering amplitude have well-defined limits and that the Coulomb screening affects only the 
Coulomb distorted waves in the initial and final states. We show how to deal with all the integrals containing the 
Coulomb off-shell T-matrix. The obtained equations are suitable for calculations and the results will be presented in 
the following up papers. 

The advantage of the developed approach is that, owe to the explicit including of the Coulomb interaction, it can 
be applied for analysis of the deuteron stripping on heavier nuclei. Such reactions provide a unique tool to study 
{n, 7) processes on exotic nuclei, which are important for nuclear astrophysics and applied physics. 

The plan of this paper is as follows. In Sect. [XX] we first derive the matrix three-body equations for the (A+2)- 
body operators projected onto target states. In Sect. IHII bv choosing (quasi-)separable ansaetze for the multichannel 
potentials these three-body equations are converted into effective-two body equations in the usual manner. In Sect. 
IIVI the AGS equations in the Coulomb distorted wave representation are derived. In Sect. [V]we present the final 
expressions for the modified AGS equations after angular momentum decomposition. In Appendices A-D the partial 
Coulomb scattering wave functions, regularized matrix elements sandwiched by the Coulomb distorted waves, the 
Coulomb off-shell scattering amplitude, the Coulomb modified form factors and the pole singularity of the exchange 
Coulomb triangular diagram are considered. Throughout the paper the consideration is done in the center of mass of 
the three-body system a + f3 + 7, that is the sum of momenta of all three particles is always zero. We use the system 
of units in which h = c = 1. 



II. REDUCTION OF THE (A+2)-PARTICLE TO A THREE-PARTICLE PROBLEM 

In this section we consider the reduction of the system consisting of two nucleons (denoted as particles 1 (proton) and 
3 (neutron)) and a nucleus consisting of A nucleons (particle 2) to the three-body system. Presentation here extends 
the formalism presented in |8j by including the proton-target Coulomb interaction and removing some inconsistencies 
in 1. 

The Hamiltonian is given as 

H = H int + H +V, (1) 
H = K 2 a /2(i a + Ql/2M a , (2) 

V=Vi+V 2 + V 3 , (3) 

where Hi nt — Ti nt + Vi n t is the internal Hamiltonian of nucleus 2; Tj„ t and Vi n t are the internal kinetic energy operator 
and internal potential of nucleus 2. Hq is the Hamiltonian of the relative motion of the non-interacting particles 1, 3 
and the center of mass of particle 2. That is, K Q is the momentum operator for the relative motion of particles /? and 
7 and fi a = mprrij/mp-y the corresponding reduced mass, m / 3 7 = mp + m 7 ; Q Q is the relative momentum operator 
for the motion of particle a and the center of mass of (/3,7) with M a = m a m^ 7 /(m Q + mp + m 7 ); m v denotes the 
mass of particle v. The potentials Vi and V3 describe the interaction of the nucleons 3 and 1, respectively, with 
each of the constituents of nucleus 2, and Vi is the internucleon potential. Potential V3 — V 3 S + Vf, where and 



4 



V§ are the short-range and the Coulomb part of the proton-target interaction, respectively. The Coulomb potential 
depends only on the distance between the proton (particle 1) and the center of mass of nucleus 2. For simplicity, in 
this section we disregard the Coulomb interaction, which will be explicitly included in the next sections. 

Consider the case that nucleus 2 can exist in several internal states p (p = 1,2 ... ,N), assumed to be orthogonal, 
with wave functions \(p% > and energies e p > 0, that is, 

ffinM (4) 

The notation is such that p = 1 corresponds to the ground state with e 1 = 0. The index p is supposed to contain 
the complete specification of the internal state, in particular also its spin, isospin etc. In concrete calculations it is 
not necessary to limit oneself to genuine bound states; also so-called quasi-states simulating the contribution from the 
continuous spectrum of the internal Hamiltonian of the nucleus might be included among the {\<f2 >}• 

To reduce the (A+2)-particle problem to the much simpler three-body problem we project all operators acting in 
the (A+2)-particle space onto the three-particle space. In this way they become N x N matrix operators: 

H = [H p °] = [< rf\H\<f% >], (5) 
H = [H p a ] = [< rtmtf >} = [S pa (Kl/2p a + Q 2 j2M a )}, (6) 

Ya = \y£ a ] = [<<rt\V a \rf>], (7) 

[Spa Q a ] = [< ^IQcJ 1 ^ >] an d [^pjKo] = [< <fi2\K a \tp% >]. The resolvent matrices corresponding to the restricted 
full and free Hamiltonian matrices are 

Q(z) = (z-H)-\ (8) 

g (z) = (z-R )-\ (9) 

All the operators acting in the A + 2 space and projected onto the three-body space are underlined, while their 
matrix elements, which have upper indices characterizing the excited states of nucleus A, arc not. Note that since the 
interaction between the nucleons 1 and 3 does not depend on the internal state of nucleus 2, the potential matrix V2 
is diagonal: [V^* 7 ] = [S pa V2], V2 being a scalar function. 

Next we introduce the channel Hamiltonian matrix for channel a: 

t a = [HS'] = [(Kl/(2^ a ) + V a ) p,T ]. (10) 

The plane wave |q Q > is eigenfunction of the operator Q a to the eigenvalue q Q and 

Q 2 J(2M a )\q a >=q 2 J(2M a )\q a >. (11) 

Let us introduce the bound state of particles (3 = 2 and 7 with the quantum numbers collectively denoted by n a 
with the wave function ip a n a , a^2. It satisfies 

(E a n a — H a — Hint) fan a = 0, (12) 

H a = K^/(2/i a ) + V a and E ana < is the binding energy of the bound state {a, n a }. Multiplying this equation 
from the left by the bound state wave function ip p 2 of nucleus 2 in the excited state p and inserting ^ \<p\ >< (p 2 | we 

cr 

obtain the coupled equations 

E^riCn. >=K n JKn a > (13) 
a 

for the overlap functions of the bound state wave function (p a n a an( i the bound state wave function of nucleus 2 

\Kn a > = < <P2\<P*n a >, P = 1,...,N. (14) 

A priori, we have infinite set of the coupled overlap functions, but here we restrict the number of the excited states 
of nucleus 2, constraining, correspondingly, the number of the coupled overlap functions by N. Also 

E = q 2 J(2M a )+E p na +e» = q 2 J(2M a )+E ana , a ± 2, 

E = q 2 J(2M a )+E ana +e p , a = 2. (15) 
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is the total energy of the three-body system a + /3 + 7 with the pair a in the bound state n a and 

EL a = E ana -e"<0. (16) 

Assume that /3 = 2 and 7 is a nucleon. Then E£ n is the binding energy for the decay of the bound state (/37)n a — > 
/3 P + 7, where f3 p = 2 P is the nucleus 2 being in the excited state p. Here we denote E anct = E^ n , a^2. For a = 2 
E a n a is the deuteron binding energy. 

Note that from normalization < ip an ' a |<v9an Q >= <5n' Q n Q we get 

N 

< <fan' a \<Pan a >= X! < ^n'J^«o >= ( 17 ) 

(7=1 

For the two- nucleon subsystem a = 2, Eq. (|13p reduces to the familiar eigenvalue equation since both nucleons are 
structureless. 

We define a matrix of two-body T-operators t a (z) in the three-body space for subsystem a 7^ 3 via the Lippmann- 
Schwinger equation 

t a (z)=Y a + Y a G (z)t_ a (z)- (18) 
This amplitude is obtained from the standard nucleon-nucleus scattering amplitude 

t a (z) = V a + V a g {z)t a (z) (19) 
by projecting it onto the target A bound states, where 

Qa{z) = ; r 1 h (20) 

Z — £Lq — Hint 

is the free Green function containing the internal Hamiltonian Hi nt of nucleus 2. 

The elements t p °{z) are related to the corresponding ones of the matrix of two-particle T-operators t Q read in the 
two-particle space as 

= / (^fjlP« > i p °{z-pl/2M a ) < p a |. (21) 
Equation (|19p is equivalent to the following coupled system of Lippmann-Schwinger equations for the operators t p °{z): 

N 1 

fr{ z a -e T - K£/(2^ Q ) 

where the energy shift accounts for the different reaction thresholds due to the excitation of the nucleus. On the 
energy shell, p a — q a , z — E and 



Also 



with 



z a = z-ql/(2M a ). (23) 
*a(2a) = [S P Mh)] (24) 



Uh) = V 2 + V2 — 2 t 2 (z 2 ) (25) 

^2 - K|/(2^ 2 ) 

being the purely elastic nucleon-nucleon T-opcrator. 

After we have defined all the necessary ingredients we can introduce the transition operators satisfying the AGS 
equations. In this paper we use the modified transition operators. To explain this modifications we consider the 
standard transition operators @j: 

U 0a (z) = -S 0a (H + T int -z)+V'-V a - V int - S 0a (Vp + Vint) + V, 3 gV a . (26) 
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Here, Sp a — 1 — 5p a is the anti-Kronecker symbol. 

V = V a + Vp + V 1 +V inU (27) 

where Vint is the interaction potential of nucleons of target A, V a = V — V a — Vint- Let us first consider the 
nondiagonal transition operators (/3 a): 

Up a (z) = -(H + T lnt - z )+v'-V a - V int -Vp- V int + V f3 QV a = -(H Q + H in t — z) + V — V a — Vp + VpQV a , 

(28) 

where the full resolvent G(z) = (z — H)~ l . Similarly for the diagonal transition (J3 = a) 

U aa = V a + V a GV a . (29) 

Note that the on-shell reaction amplitude is given by the matrix element < $^ n(3 \Up a (z)\<S> ana >, where the channel 
wave function |$ a „ Q >= fane |0a >j lOa > i s the plane wave in the initial channel a. Because (ff + H int + V a — 
z) &an a = 0, the matrix element < $p np \Up a (z)\$ ana >=< &p nfi \Up a \z)\$ ana >, where U^(z) is the standard 

transition operator U^~\z) = V p + V p GV a . 

Thus we can rewrite the standard transition operator in the form, in which Hi nt can be combined with Hq as a 
"modified free motion" Hamiltonian Hq + Hint- After that we project the transition operator Up a onto the eigenfunc- 
tions of Hi n t, what allows us to eliminate from the consideration Hint replacing it by the corresponding excitation 
energy of the target e p . The modified transition operators U_p a (z) also satisfy the AGS equations |3| (the Coulomb 
interaction between particles 1 and 2 is, for the moment, disregarded) : 



U# a {z) = SpaG^iz) + £M^ 7 (20go(*)M Z )' ( 3 °) 
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It is apparent that the general form of ([30|l coincides with that for three point particles, the only difference being 
that in the present case all operators are now N x N matrix operators. The physical amplitude for the transition 
from an incoming state in channel a with relative momentum q Q , with the bound state of particles /3 and 7 being 
characterized by quantum numbers n a , to an outgoing state in channel /3 characterized by relative momentum qp and 
bound-state quantum numbers np, irrespective of the internal excitation state of nucleus 2 and taking into account all 
intermediate-state excitations and de-excitations of the nucleus due to scattering with each of the nucleons, is defined 
as matrix element of U_ a p between the channel states $ „ oiqo : 

E + iQ) =< <$>p n ^^\Up a {E + iQ)\<5> ana ^ a > (31) 

with E being the total energy of the system. The introduced channel wave functions |3>cm ;q a > have the form of a 
column matrix, 

>= [K„ > |q Q >] = { . > ■ (32) 

Equation (131 [) represents the sum of all contributions from all excitation states of the nucleus in the incoming and the 
outgoing state which are allowed by the on-shell condition 

E = q 2 a /2M a + EP na +6? = q 2 a /2M a + E aHa 

= q'p 2 /2Mp + E$ np +e° = q f/2Mp + Ep np . (33) 

The component 

X 7n„ a nM^E + iQ) =< q'p\ < ^ | (E + «0) | C» Q > |q« >> (34) 

where \Up^\ = [< ¥?f |^/3q|<^2 >L describes the transition from an incoming a-channel configuration {q Q ,n Q }, with 
the nucleus being in excitation state p and satisfying the on-shell constraint, to an outgoing channel /3-channel 
configuration characterized by {q'^,np} and internal excitation state a, taking into account all intermediate-state 
excitations and de-excitations of the nucleus due to scattering with each of the nucleons. q Q and q^ are the on-shell 
relative momenta of particles in the initial channel a and final channel /3. 
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III. SEPARABLE MULTICHANNEL POTENTIALS AND EFFECTIVE-TWO BODY EQUATIONS 



As is well known, the solution of the AGS equations (|30|) is greatly simplified if the subsystem transition operators 
t a are represented in separable form. For the nucleon-nucleus subsystem 2) we assume the interaction V£ a 
leading from a two-body state with internal excitation state p of the nucleus, to a two-body state with internal state 
a of the latter, to be described by a (quasi-separable) multichannel potential of the form 

A a 

V n[n a = E ^ X an' a t' a > K^-n' a t' a n a t a < Xan a tJ ( 35 ) 
t' a t a 

for p, a = 1, 2, . . . N. Here, n a (n' a ) denote the quantum numbers of the two-body state in the pair a before (after) the 
interaction, t a (t' a ) is the number of the separable expansion term before (after) interaction with the total number of 
the expansion terms A a . A priori, A a > N a , where N a is the number of the bound states in the pair a. Terms with 
t a > N a represent auxiliary terms, which don't correspond to bound states but are needed to improve the accuracy. 

The form factor for the state with quantum numbers n a , term number of separable expansion t a and inter- 
nal excitation a is denoted by \Xan a t a >■ The coupling matrix [X^ n , t , n t ] is chosen to be symmetric, that is, 
Kx-ri v n a t a ~ Kx P n a t a n' v m or der t° ensure the hermiticity of the potential. Since, as mentioned above, the nucleon- 
nucleon interaction does not depend on the internal state of the spectator nucleus and is therefore diagonal over the 
upper-scripts, we have for the matrix elements of the coupling matrix 

^2°n' a t' a n a t a = tifxr^2;n' a t' a n a t a - (36) 

The solution of the Lippmann-Schwinger equation (|22l) for the above potential matrix (|35|) has the form 

A a 

C(^) = E \X*<t' a > K- K t >a t a (2«) < Xan a J, (37) 
t' a t a 

where the A^, t , nata {z a ) satisfy the coupled equations 

N A a ^ 

^a;n' V n a t a (^<*) = ^a;n' t' n a t a ^°' :n ' a t 'a n 'a t 'a ^ ^««!J*!i I a _ c t _ JS2 /o,, IXan"t£ > ^a;<t>„i„(^a) 

t — 1 n" t" 

(38) 

for a 2, with a similar equation for a — 2 which is, however, uncoupled with respect to the indices characterizing 
the internal state of the nucleus. 

Now we will derive the expression for the overlap function | (f> p ana > in the separable representation. To this end we 
can rewrite Eq. (|13[) in the more convenient form: 

(Kn a \<n a >= E VST\<f>' > • (39) 



2 Pa 



From this equation taking into account Eq. (|35l) for the multichannel potential we get for the overlap function with 
quantum numbers n a of the pair a^2 and with the nucleus being in excited state p 



an a 



where 

JY 



■^KT E C Lt a \Xan ata >, (40) 



= E E E K-n aWa < Xln> a t> a \4>ln> a > ■ (41) 



From normalization condition (|17[) we get 



N Aa r 1 i2 

E E O a <A <xL at j[ . p _ K , ] lxL Qf „>=l- (42) 
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The denominator in this equation is nonsingular because E p na < 0. Eq. (1401) is extremely important for solution 
of the Faddeev equations for deuteron stripping in the separable representation. The overlap function is not an 
eigenfunction \4> p IJla > of any Hermitian Hamiltonian and, hence, not normalized to unity. Its square of the norm is 
the spectroscopic factor of the configuration f5 p + 7 in the bound state a„ a : 

< <» a Kn« >=%"7)n a - ( 43 ) 

In the three-body model, in which the antisymmetrization between the nucleon 7 and nucleons of nucleus /3 = A 

N 

is neglected, the sum rule is S^p P ^ n =1. If we take this antisymmetrization into account, for what we need 

P =i 

N 

to go beyond of the three-body approach, in the isospin formalism the sum rule is ^ Srpp 7 -\ n = A + 1. For 

P =i 

N 

antisymmetrization separately with respect to protons and neutrons the sum rules is ^2 S^ P ^ n = Na + 1 for 7 = 3 

N 

and ^ S{pt"fin = %A + 1 f° r 7 ~ 1j where Na {Za) is the number of neutrons (protons) in nucleus A. To take into 
P =i 

account the antisymmetrization we can include the antisymmetrization factor into the coefficients c^ q4q . 

Because we use the three-body approach we suggest the following procedure. For a / 3 the overlap function v?an a 
has a pole at k\ = — 2 // Q E p n . Its residue in the pole can be expressed in terms of the asymptotic normalization 
coefficient (ANC) for the virtual decay a na — > f3 p + 7. The neutron ANCs can be determined from the analysis of 
the experimental data. The overall normalization due to the antisymmetrization can be included into the coefficients 
c p na . For a = 3 the singularity at k 2 a — ~1[i a EP, n is the branching point but it is a pole for the two-body p + A 
scattering T-matrix and its residue is expressed in terms of the proton ANC. For a = 2 the overlap function is the 
deuteron bound state wave function, which is decoupled from the indices characterizing the excitation of nucleus A. 

Now we proceed to the generalized AGS equations. Let — {n a , t a } collectively denotes the quantum numbers of 
the state of the pair a and the number of the separable expansion term. Allowance for the Coulomb proton-nucleus 
interaction leads for the physical transition amplitudes, which are the matrix elements of the operators W^S: 

^£ f ,aC^>1'i z ) = < ^1 < 9?cJ5o' WW^IflSc. > |q Q >, tp<Np, t a <N a , (44) 

where for a = 3 is the Coulomb-modified form factor, see Eq. (|84|) below. For a^3 9 a c a ~ ^a( a - The total 

physical amplitude for transition from channel a with bound state quantum numbers n a and internal state p of the 
nucleus A to channel (3 with a configuration characterized by na and a is determined by 

N fi N a 

xz P , a nM^^;z) = E E m, a cM^ z )- ( 45 ) 

Note that 

ew< f ,, a cJ*) = Sf>a*%c> aCa (z), (46) 
where £/9 — C a - The amplitudes A^ a ^ satisfy the off-shell effective- two body equations 

Ay N f dp 

X !c<3,«Co( P /3' P " ;2: ) = ^P Z ^,aCo(P/3'P«! Z ) + ZJ X, X, X, / J5~~)3 Z P(( l ,~fC,(P'p'>Pj'-> Z ) 

7 n^n^ ti r t 1 = l r=l 

x ^C^-pV^^I^cAp^P^z)- (47) 

Here and p a are the off-shell momenta. The off-shell amplitudes are needed later on when deriving the AGS 
equations in the Coulomb distorted wave representation. 

For only two charged particles (with the charges of the same sign), particles 1 and 2 in the case under consideration, 
the effective potentials are given by [13] 

zrc P ,«M>P«'> z ) = s°e$Z%*<Pp\<tf<^ >\p a >- (48) 

Here the Coulomb Green function operator is 

= Bo "Y?)" 1 (49) 
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and \Q% aa } = [< <p£|0f K >]> Ws" 9 ] = [tvpVf]- The effective potentials are diagonal in the upper indices since 
the internal state of nucleus 2 cannot be changed by the one-step particle transfer process. Equations (f4"?| allow 
one to take into account the possibility of the excitation of nucleus through the interaction with the nucleons 1 and 
3, respectively, as well as the contribution of the rescattering of the excited nucleus to the three-body dynamics. 
Because our Coulomb V3 interaction potential depends only the distance between particle 1 (proton) and the center 
of mass of nucleus 2, this Coulomb interaction cannot excite nucleus and [V 3 ap ] is diagonal matrix in the upper 
indices indicating the target excitation. 

Further insight is gained on expanding Eq. (|47[) in a Neumann series. If in that expansion one neglects the con- 
tribution from all those terms in which the excited nuclear states appear explicitly, one obtains equations which are 
of the same form as the usual AGS equations. The only difference between the latter and the AGS equations for 
three point particles is that now the two-particle amplitudes describing the elastic scattering of each of the nucleons 
off nucleus 2 take into account the multistep excitation and subsequent de-excitation of the nucleus. Therefore, the 
corresponding pair potentials must be expected to be energy-dependent, non-local, and complex. In practical appli- 
cations these potentials could be approximated by optical-like potentials but fitted to the nucleon-nucleus scattering 
and bound-state data. 

IV. THREE-BODY EQUATIONS IN THE COULOMB DISTORTED WAVE REPRESENTATION 

In this section we obtain the AGS equations in the Coulomb distorted wave representation. To this end we start 
from the matrix AGS equations (|47p. a ^ {p'piMai z) is the reaction amplitude describing the transition from 

the initial channel a + (/3 7) f Q and nucleus A in the internal state p, whichever particle it is in the initial channel, 
to the final channel /3 + (a 7)^ with particle A in the internal state a. If, for example, a = 2, the above transition 
is a p + (f3 7) £ a — > /3 + (a CT 7)^. Because we use separable potentials, the effective potentials a £ Q (P^ Qa; z ) 
are sandwiched by the separable form factors XaC an< ^ X^a ' ^ n ^ e nex * section the full AGS equations with 
explicit indication of spins and angular momenta will be presented. It is worth mentioning that we assume that all 
the effective potentials are diagonal over the upper-scripts p and a, that is the effective potential doesn't change 
the internal structure of nucleus A. This change can occur only in the separable potential matrix A^^, * . This 
assumption is justified because the Coulomb p — A interaction depends only on the distance between the proton and 
the center of mass of target A. 

Our goal is to obtain the modified Faddeev equations in the AGS form where the Coulomb rescattering of the 
particles in the initial and final states is explicitly taken into account by sandwiching the transition amplitudes and 
effective potentials by the corresponding Coulomb scattering wave functions. That is what we call the Coulomb 
distorted wave representation. To obtain this representation we follow the strategy outlined in [16j. First, we need to 
rewrite the matrix AGS equations in a conventional Lippmann-Schwinger form X = Z + Z Gq X , where Go is a two- 
body free Green function. After that we add and subtract the Coulomb potential between the particles in the initial 
and final channels introducing Z = Z' + U c , where Z' = Z — U c . After that we will apply the two-potential equation 
leading to < p c \X\p c >=< p c \Z'\p c > + < p c \Z' G c X\p c >, where G c is the two-body Coulomb Green function 
describing the propagation of the particles in the intermediate state and p c is the Coulomb distorted wave. Using 
the spectral decomposition of the Coulomb Green function G c we immediately arrive at the desired AGS equations 
in the Coulomb distorted wave representation, where the reaction amplitudes and the effective potentials in all the 
channels a/3 are sandwiched by the Coulomb distorted waves: < p c \X\p c >=< p c \Z'\p c > + << p c \Z'\p c > 
Go < p c \X\p c >>. We have only quite schematically described our strategy. Its practical implementation has many 
peculiarities which should be overcome. Now we proceed to the practical implementation of the outlined strategy. 

First, to rewrite the AGS equation in the Lippmann-Schwinger form, we introduce new effective two-body amplitudes 
and effective potentials [l6j : 



X 




-1 



(50) 



and 



Z 




(51) 



with 



where C,p = Co 



z a = z- p 2 a /(2 M a ). For t a < N a 



p«) a; 



Pa) gO; a ( a {Pa; Z 



(52) 



(53) 
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<7o ; .CAP*,*- E ma ) = z „z ana - P l,2M a (M) 

is the two-body free Green function, z — E a n a — Pa/2 M a — 9a/2 M a — p a /2 M a + i 0. 

Introducing new reaction amplitudes and effective potentials we are able to rewrite half-off-shell Eqs. (|4"7| in the 
matrix Lippmann-Schwinger form: 

^ N f dp 1 

7 n T t T = l r=l" y V ; Z ^7 / V z iW 77 

x ^,ae-CP7,q«;*) + E E EE/ 72fT3^c^7c;(p^P7;-)A^c»(--P7/( 2 ^7) 



x X 7 ^„ QCa (p 7 ,q Q ;z). (55) 

Note that the two-body free Green function appears only in the terms describing the physical 
bound state — > bound state transitions. The new effective potentials Z are not diagonal over the upper indices 
because they contain now Q ap . 

The terms with t' , t" > N 7 + 1 are auxiliary terms included to improve the accuracy. Corresponding auxiliary 

amplitudes X^,, (p 7 , q Q ; z) don't describe any physical transition bound state — ¥ bound state. If we include these 
auxiliary amplitudes, we need to supplement AGS equations by the equations for these amplitudes. 

In Q the Coulomb p — A interaction was neglected but here we explicitly include it. Allowance for this Coulomb 
interaction leads to the appearance of the Coulomb-modified separable form factors for the s yste m (1 2) and the off- 
shell 1 + 2 Coulomb scattering amplitude in the four-ray vertex in the triangular diagrams pl ll7l [l8j . The triangular 
diagrams describing the d + A and p + (n A) elastic scattering with the four-ray vertex, owe to the presence of the 
off-shell p+A Coulomb scattering amplitude, contain the Coulomb forward singularity A~ 2 in the transfer momentum 
plane at A 2 = 0. Coincidence of this singularity with the pole singularity of the two-body Green function leads to a 
noncompact singularity of the generalized Faddeev equations written in the AGS form [j, 0, [lj| . To eliminate it in 
4, 16] the channel Coulomb potential was added to the diagonal effective potentials and subtracted: 

^SCa.aCa = <W<L3<$C/3Cc, <W U£ + Z^ a( . al t p or t a <N a , (56) 

and 

, a Ca = % ,a ia , tp>Np + l Ort a >N a + l. (57) 

Note that the nondiagonal effective potentials are not affected by the Coulomb channel potential, that is nondiagonal 
potential Zp£ a Q coincides with the original effective potential Z^ a( . , while its diagonal part in the channel 
a 7^ 3 is Z^ a( * — Jj£ for t a < N a ; is the channel Coulomb potential describing the interaction of particle a 
and the system (/3j) with its charge concentrated in its center of mass : 

V C {pa) = (5g) 
Pa 

Z a e is the charge of particle a and e is the charge of the system (^7), which is equal to Z$ e or Z 1 e depending 
on which particle of that system has nonzero charge; p a is the distance between a and the center of mass of the 
system (/3j). The explicit expression for the effective potentials will be given in the next section. 

After adding and subtracting the channel Coulomb potentials, according to the outlined strategy, we can apply 
the two-potential theorem, which allows us to rewrite Eqs. (|55|) in the form, in which the reaction amplitudes and 
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potentials are sandwiched by the Coulomb scattering wave functions (off- and on-shell): 

N -> N r a 

^Pd^tJPp > q Q > z J - ^c/s^CoAP/s ,q Q + 2^ 2^ 2^ 2^ y ( 27r )3 ^o^c^P/? P7 ^ z > 



7 n 7 i T — 1 r— 1 ' 



1 x sc 1 r,'"'- 1 



2 



£ 7ni -p 7 V(2M 7 ) 



2/(0 M\ " 1 7Ct,"Co^7 '1" 



+ E E E E / 7^f^( p ; c( ^ p ^^^^^^ 

tp < N fi , ta < N a , (59) 

fSCap / / _\ _ y'SCcrp ( I C(+l \ 

A )3( fl ,aC«lP/!' > z ) ~ Z /3 Cl 3 ,« Cc ( P /3 ' q " ' ^ 

W-r JV , 

, V V V / Pl 7SCar ( I „C(-). n i j>SCrp ( C(-) C(+) . \ 

+ 2^ 2-, 2-, 2_, / f 2 7r1 3 Z /3C^7C^P/3'P7 _ n 2/rtMl A 7C„oC^ P 7 ' <la . *J 

7 n-, fc y =lT=l J V ; Z ^yrh ^7 l\* m l) 

+E E E E / w&zr^M'ih-.*) ^cM-vV^^^lZtS^ <£ (+) ;-)> 

^ >JV^ + 1, t a <N a . (60) 

These equations are the desired Faddeev equations in the AGS form in the Coulomb distorted wave representation. 
They generalize Eqs. (5.35) [l(| by taking into account the target excitation. Eqs ([59)) determine the components 
of the reaction amplitude matrix corresponding to bound state — » bound state transitions. From these components 
one can calculate the observable cross sections. However, to determine these physical matrix elements we need to 
know also the auxiliary amplitudes X^, T ^ (p 7 , q« ; z) corresponding to transition from the initial bound states 

(t a < N a ) to the quasiparticle states (t' > iV 7 + 1). That is why Eqs (|59"T) are supplemented by Eqs (f60|) . which 
determine these auxiliary components. We remind that the auxiliary components appear as the result of the separable 
expansion of the nuclear potential not only over the bound states but also over the quasiparticle states not representing 
any physical bound states (tp > Np + 1). 

In these equations the reaction amplitudes are 

*IST£c>*«£ (+) ;*) =< y/>l%.«c.(*)l<. (+) >> ( 61 ) 



where for tg < Nr. t a < N a yg = p R ( ' = ib , ) denotes the off-shell Coulomb scattering wave function in the 



exit channel (3 and for tp > Np + 1, t a < N a yp = p^; q Q — Vqi stands for the on-shell Coulomb scattering 
wave functions in the initial channel a. Also for the amplitudes under the integral 



ySCtp 



X 



7 cw>^ (+) ;*) = < y-r\x T ,l, a cM^ +) >. ( 62 ) 

c(-) Of— ) 

where for t 7 < A 7 , t a < N a y 7 = p 7 — ip p } f is the on-shell Coulomb scattering wave function in the interme- 
diate channel 7 and for t 1 > N 1 + 1, t a < N a y 7 = p 7 . 
The inhomogeneous terms are 

^i c ; p c „(py ( ^q« (+) ;-) =< ^:i\Kz^M^ +) >> h < *« < (63) 

%TaC.(Pfl' < ^ (+) 5 Z ) =< Pp-|%,aCc(*)l<i +) >' ^ > ^ + 1, *a < N a . (64) 

The potentials in the integrand are 

z; s ^i(pf { -\ P f-hz) =< ^l\z'pZ, 7C M€t } >' ^ N ^ *7 < (65) 
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^l^^^^^^^l^^cM^^ h<Np, t 7 >JV 7 + l, (66) 

^fS"tr(PP'P7 (-) ;*) =< P^W^IC" >< */» > ^ + 1, ty < JV 7J (67) 

^, 7 ^(P^P7^)=<P^|^,7^W|P7>> ^>^ + l, t 7 > AT 7 + 1. (68) 

In each channel, for which two-body state is not a bound state, the Coulomb scattering wave function (off-shell or 
on-shell) should be replaced by the corresponding plane wave. 
The off-shell Coulomb scattering wave function is given by 

KSi > = t 1 + Go(«£/(2M a )] TC(q 2 a /(2M a ))\ Pa >, p a + q a , (69) 

G (q 2 a /(2 M a )) = 2/r9 „, 02/^9 «f^-n (70) 
qil (2 M Q ) - Q£/(2 M a ) + * 

is the free Green function describing the propagation of the system of the noninteracting particle a and the system 
(/?7), Q^/(2M a ) is the kinetic energy operator of their relative motion. Also T^(p' a , p a ; ^/(2M„)), p' a , p a ^ q a , 
is the two-body off-shell Coulomb scattering amplitude of particle a and the center of mass of the system [j3 7) 
moving with the relative kinetic energy q^/(2M a ) and interacting via the Coulomb potential U^(p a ). 

The on-shell Coulomb scattering wave function can be obtained from it by taking the limit p a — > q a . In this limit 
we get [H 

lim <^7l\n(p a ,q a )=<^->\, (71) 
Pa— >q a r 1 ^ 

where 

fi(p«, q a ) = ^ W2 - iVa )}- 1 ( Pa+qa . n Y n \ (72) 

r/ct = Z a Zp 1 e 2 M a /q a is the Coulomb parameter in the channel a, which characterizes the Coulomb interaction 
between particle a and the system /3 + 7 moving with the relative momentum q a , Zp~ = Z p + Z„. 

We return now to Eqs (f5T)l) . These are not integral equations yet and we will address this point. Let us consider the 
physical amplitudes (describing the transition bound state — > bound state) . On the left-hand side of Eq. (|59| we have 

the amplitude Xp^a^Jj? 

;z) given by Eq. (JBTJ, while in the integrand we have different half-off-shell 

amplitude A^£^ c Jp ; z) given by Eq. (JH3) - Both amplitudes are half-off-shell but the off-shell effects are 

treated differently in both amplitudes. In X^"^ (p^g ) Qc* ^ ; z) in the bra state we have the off-shell Coulomb 

scattering wave function, while in X^ T ^ (p 7 qa ; z ) in the bra state we have the on-shell Coulomb scattering 
wave function but with momentum p 7 , which is the integration variable and, hence, p 7 7^ g 7 , where q 7 is the on-shell 
momentum in the channel 7. Hence the amplitudes in the left-hand side and in the integrand of Eqs. (|59p are not 
the same functions and these equations cannot be solved as integral equations. 

Let us take the on-shell limit p'p — > q'p in Eqs (I5T)1) . Taking into account the on-shell limit of the off-shell scattering 
wave function, see Eq. (|7T1) . we get for the reaction amplitude and effective potentials: 

, lim , xllZcSPf^^hz) = W fi , ^T'xIlZtJ^-^^hz), (73) 

; im ^ c j Ca (p'^-\cf a M-,z) = w» ^)]- i ^c.(i 0( " ) »^ (+) !*)> ( 74 ) 

lim, Z' p %^^-\ V °^;z) = W„, <t p )]-lZ™°^ffi-\ V y-\z) (75) 



and 



lim, Z s ^;^(p'^-\p Y ,z) = W p , q'pT 1 Zfly c ^-\p Y ,z). (76) 
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Thus taking limit p'p —> q'p in Eqs ((59)) and multiplying them by [Q(p'p, q'p)] we obtain the on-shell limit 

jrSCvp ('C(-) C(+). \ _ y'SC ap ('C(-) _C(+).-A 

N '< N r A i 

+E E E E / ^.(^^ P^)A;: (l(; (,-^/(2M 7 p^: (jP7 , q ^;4 

tp < Np, t a < N a . (77) 
The on-shell limit of Eqs. (|60[) is straightforward because the bra state is the plane wave and we get 

\" \" \" \" / °-P7 ySCVr / / C(-).„\ i jrSCrp ( n C(-) „C(+). y \ 

7 n T t 7 =l r=l ,/ y ' Z ^7 /l ZJ "7^ 

A 7 iv 

+E E E E 

7 n. 

tfi>Np + 1, t Q < JV a , (78) 



E EE/ tI^^c^p^^ q£ (+) ;^ 



where £ 7 = {n 7 , i 7 }, £ 7 = { ra 7 i C}- Evidently that Eqs (J77J) are not integral equations because on the left-hand 
side we have the on-shell transition amplitudes while in the integrand they are half-off-shell. To obtain an integral 
equations from (1771) we use its off-shell extension, which differs from the one used in (|59[) : 

fSCap ('C(-) Cf(+) . _\ _ 7 SO (T P ('C(-) „(?(+). _\ 
A 0S/3, a ( a (Pp ' i z ) - £ Ptp,at a \Pp ' q « ' I 

N y N 

\" \^ \^ \^ / d P7 d-'SCar ('C(-) ± jrSCrp /_C(-) „C(+).^ 



(2*)» "^-P ' '.-I^-^/(2M 7 ) 

^7 N 

EE EE/ ^^S^(P? (-) . P,^)^c;(^))^c;ZcJP- q 2 



^ J (2tt) s 

tp < Np, ta < N a , (79) 



ySCap ,1 . _\ _ 7 SC op ( , C{+) \ 

A pCp,*to: [ Vpi Qa. ,Z) - ^Pi^aCc^P'^c > Z > 



W V V [ AVl z'SCar ( I C(-) ) 1 j>SCrp ( C(-) C?(+) . 

Ve^^i (2^) 3 ' Z) z-£ 7 „ 7 - P7 V(2M 7 ) A 7C 7 ,ocJP7 - q Q - 



t_ = l r=l ' 



2 



^7 iV . , 

+ E E EE / Tfe^c,^^;*)^^^))^,^^ q ^ (+) ^)^ 

*p>^ + l,t tt <JV a> (80) 

where z 1 — z -p 7 /(2M 7 ). Eqs (|79p are half-off-shell because in the exit channel /? the on-shell Coulomb scattering 

wave function < | =< ip^i \ is present with momentum p'^ ^ q'p. It means that in Eqs ([79")) the off-shell 

behavior in the exit channel differs from the one in Eqs ([59)) , where the off-shell Coulomb scattering wave function is 
used rather than the on-shell one but with the off-shell momentum. But the on-shell limit of (|T9")) coincides with the 
renormalized on-shell limit of (|59|) . Eqs (|79|) together with (l80l) are our final equations in the Coulomb distorted wave 
representation, which will be used to calculate the deuteron stripping amplitudes and cross sections. The advantage of 
these equations is that one don't need to use Coulomb screening procedure, application of which becomes problematic 
when the charge of the target increases. 
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V. ANGULAR MOMENTUM DECOMPOSITION OF GENERALIZED AGS EQUATIONS IN THE 

COULOMB DISTORTED WAVE REPRESENTATION 



Here we present the final expressions for the modified AGS equations |79|) and (l80|) after the angular momentum 
decomposition. We follow the formalism used in [l5j . We use the following angular momentum coupling scheme for 
a given channel a: + s 7 = S a , L a + S a = J Q , s a + J Q = E Q , l a + S Q = J. Here, s Q denotes the spin of 
particle a, L a is the relative orbital angular momentum, S a is the total spin, and J Q is the total angular momentum 
of particles /? and 7; moreover, 1 Q denotes the relative orbital angular momentum of particle a and the pair (^7), 
and finally J is the total angular momentum of the three-body system. Also A a = A JaSa is the rank of the separable 
expansion in the two-body channel a with fixed J a and S a ; N a = N JaSa (7V Jq,Scv < A JaSa ) is the rank of the 
separable expansion corresponding to the bound states in the pair a for given quantum numbers J a and S a ; as before 
t a = l...A JaSa enumerates the number of the expansion term of the separable potential in channel a with given J a 
and S a . As in the previous sections, n a = {L a , S a , J a } collectively denotes the quantum numbers of the pair a, 
while £ Q = {n a , t a } = {L a , S a , J a , t a } denotes the complete set of quantum numbers characterizing the two-body 
state a, which, in addition to n a , includes the number of the separable expansion. A new variable u a — {n a , l a , S Q } 
collectively denotes all the introduced above spin-angular momentum variables in the channel a except for t a . As in 
[1] we don't use the isospin formalism to treat two nucleons. 

From Eqs jZSJ we get 

N °° j 2 

ySC ,r up , 'C C \ _ y'SC r op , 'C C \ , / y' SC ,r ar 1 'C C \ 

^upt ff ,u a t a \Ppip> Qal a i z ) — ^uptp,u a t a \P/3lp' q a l a i z ) 2-^i 2-^i Z_* 2-^1 ~2lT^~ 13 tf> ' U '< 4-7 18 1/3 ' Pl ll ' Z 



7 U~f try—1 T—l 







I jSCJ-rp ( C C \ 

*/j < tyj, *a < N a . (81) 



Here n 7 — {S 1 , L 7 , J 7 } and n 7 = {5 7 , i 7 , J 7 }, u 7 = {n 7 , £ 7 , £ 7 } and u 7 = {n 7 , £ 7 , S 7 }. Note that in the scattering 
particles of the pair a 
Also Eq. USD) leads to 

fSCJ'ap / ' C . _ ySC Fop / ' C . \ , f d PlPj yrSCT ar ( ' C \ 

^up tfj, u a t a \Pp' 1al a ' Z ) ~ ^up tp,u a t a \P/3> Valci z ) 2-^1 L^i 2 7T 2 3 * p ' ™ 7 4-7 ^ ' 7 ' T ' 

^, ^, J 1- 1 « 



of particles of the pair a the total momentum J 7 is conserved. That is why A^ 7 "^ 7 ^ t , (z 7 ) is diagonal over J 7 . 

IV. |V 00 , o 
' dp 7 p2 

7 u T £^,=1 r—1 q 

l 



Ejn-, -p 7 2 /(2M 7 ) 

dp 7 P 7 



2^ 2 7T 2 "''H'^'i^' Pl^l^L^t-yL^t'^l) A ^t;,« ot t c «VP7' <iala'' Z )' 

tp>N{, + 1, t Q < JV a . (82) 



^■L^ty °L l v (%) sa ti snes the system of equations 



iJyS^TV f&\_ \J-,S-,TV , TW J 7 S 7 W, J- I JjSyU iJ-fS-tUlV \ 



L" t",t"' = l w=l 

7 7' 7 



^7 2 /(2M 7 ) 

(83) 



where z 7 = z — p 2 /(2 M 7 ). 

The Coulomb-modified form factor is given by 
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-l 



where 

1 

T 7 %(fc;,fc 7 ;z 7 ) = i I dxP^(x)T^(k , 1 ,k 1 ;z 7 ), (85) 

where x = k 7 • k 7 , k = k/fe. 

The physical reaction amplitude ^i^s/j/f s i' which describes the transition from the initial channel a, where the 
particle pair ((3j) is a deuteron with J a = Ji = 1, the relative orbital angular momentum of particle a and deuteron 
is l a = U and the total channel spin E a = T,i, to the final channel /3 where the particle pair (017) is in a bound state 
J/3 = Jf and the channel orbital angular momentum and channel spin lp = If and T,p = respectively, can be 
calculated from the solutions of Eqs (|81l) and (1821) as 

^/e/j^s.i^^' 1aiJ E +) = 5j e J f 5 hh fe/ss, S Ja i8i aU S^ aSi X^/^^ip^, q% la ;E + ), (86) 

S+ = S + iO. 

Now wc proceed to the effective potentials, which are the main input into the AGS equations. Our new defined 
potentials, see Eqs. (I5T1) . (|52"j) and (j53"j) . are not diagonal over upper indices. We start from generalizing Eqs (33) 
and (34) obtained in [15|, which determine the off-shell effective potentials sandwiched by the plane waves. Here we 
replace these potentials by the ones in the Coulomb distorted wave representation: 

4 

y'SCJ^op 1 'C C \ _ fr'SC J* crp{i), 'C C \ 

^up tp,u a t a \P/3lf,i Pal a i z ) — 2_*t «/3 tp,Ua *= Pal a \ z )i 

i=0 

tp < N J ? s e, t a < N J " S ", (87) 

4 



ySCJ^op 1 1 C \ _ ySC J*" op (i) 1 1 c \ 

^up tp,u a t a \Ppi Pal a i z ) — 2^1 up tf,,u a t a \P/3i Pal a i z )i 



and 



> N J v s v +1, t a < N J * Sa , 



4 

z Z?p,u a t a (P|s> P*> z) = ^ z uf, t^ul t a (p'p> p a ;z), 

i=0 

tp > N J f> s f> + 1, t a > N J " S " + 1. (89) 



Here 



y,'SCJ*op{i) ( 'c O n _ 7 d PpP0 f dPaP'a ,C fJ/w'i'^W (f Jl n »- 7 \,h C (r>"\ 
Z/ u /3 t f) ,u a t a \Pl3lpi Pal a ' z ) — — / — 2n^ — ^P'fih^l 3 ' Zj Uf j t f ,,u a t a \Pfi-> Pai z )y Pa l a \Pa)i 



tp < N J ? s e, t a < N J <* Sa , (90) 

ySC J" <jp{i) C .~\—f dPaPa_ yJ" °P (t) // // s , C //M 
ifl.tia i a VP/3' Pal a i z ) — I — 2~J^ — Zj u /1 t l3 ,u a t a \Pl3^ PaT Z )V ; Pa l c ,{Pa)y 



> iV' 7 ^ + 1, t a < 7V J ° S °, (91) 
^;,il(P^Pa^) = E^>^^)^ tp<N J " s ", t a <N J ^, (92) 

<r/S t j^,K^) = E^(^'<)<£i*J^'P-^)' ^>^ + i, i Q <tv^, (93) 
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(23) 




(13) 



FIG. 1: Pole diagram describing the neutron transfer in the reaction 1 + (23) — ► 2 + (13). 



(13) 




(23) 



FIG. 2: Pole diagram describing the neutron transfer in the reaction 2 + (13) — ¥ 1 + (23). 



and 



^ Uj3 t^ujte, (P'0) Va\Z) — A l K (p'p, Pa) R Ufj tp,Z a ta (Pp> Pa\ Z )l 



tp > N J f s f + 1, t a > N JaS " + 1. (94) 



Here 



^u/tp,u a t a {P% Pa: z ) ~ Ru/i tj,< t' a (P% Pa'i z ) ^i/ftl,Y a t a ( z a) (%a — Ean a ), 



N 



J aa (i) 



L' a t' a = l 



(95) 



N J a s a 

K P ? tlu a t a (Pp, Pa, Z ) = J2 R Ztt< t' a (P'P> ^ Z, t a ( ic ") ^ a ~ E<*n a ), 

t p > N J ^ S ^ +1, t' a <N Q , (96) 

where z a = z — p" a /{2 M a ) . The nondiagonal potential for i = describes the pole diagram corresponding to the 
neutron (particle 3) transfer, see Figs [1] and [21 the potential for i = 1(2) describes particle 1(2) transfer diagrams 
which contain one Coulomb-modified form factor, see Figs [3] and 2] The pole neutron transfer diagrams describing 
the inverse processes are shown in Figs [5] and [5] The effective potential for i — 3 describes the elastic or inelastic 
scattering triangle diagram, see Figs [7] and El and the effective potential for i = 4 describes the exchange triangular 
diagram leading to the rearrangement in the channels a^3, see Figs l9l and 1101 

Both triangular diagrams contain the 1 + 2 off-shell Coulomb scattering amplitude in the four-ray vertex. The 
amplitude of the elastic scattering triangular diagram contains a strong forward Coulomb singularity generated by 
the off-shell Coulomb scattering amplitude, which is compensated by subtracting a corresponding Born Coulomb 
scattering term in the channels a / 3. For the contributions with i = 0, 1, 2 ft = £ is a single index while for 
i = 3, 4 it is a multi-index k — (Ci, £2, /)• Explicit equations for A l K and Rulu, u a t a were given m [HI- Equations 
for A l K can be taken from this paper with only two minor modifications: q'p should be replaced by p'^ and q a by p a . 
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(12) 2 (12) 2 




3 (13) 3 (13) 

FIG. 3: The pole diagram describing the proton transfer in the reaction 3 + (12) — > 2 + (13), which contains the Coulomb- 
modified form factor (12) — ► 1 + 2. Because this form factor consists of two terms, see Eq. (|84p , the diagram is also represented 
by the sum of two diagrams. In the first diagram in the vertex form factor (12) —¥ 1 + 2 no Coulomb interaction is included, 
but it is included in the second diagram. Bubble shows the off-shell Coulomb scattering amplitude of proton 1 and nucleus 2. 



(12) 1 (12) 1 




3 (23) 3 (23) 

FIG. 4: The pole diagram describing the nucleus transfer in the reaction 3 + (12) — > 1 + (23), which contains the Coulomb- 
modified form factor (12) — > 1 + 2. Because this form factor consists of two terms, the diagram is also represented by the sum 
of two diagrams. In the first diagram in the vertex form factor (12) —5- 1 + 2 no Coulomb interaction is included, but it is 
included in the second diagram. Bubble shows the off-shell Coulomb scattering amplitude of proton 1 and nucleus 2. 



Equations for K should be modified by including the target excitation indices. In particular, Eq. (Al) [l5[ for 
i = 0, 1, 2 takes the form 

E aa ^ c ( n ' n •z)-R" r ® C M o -z)-7> a c°*c° - f Jt P. (r) fc ° ^ L %* (fc ^ 9 t ^ ,• - n 1 2 

R Cp<a {p p ,p ai z)-R Ci3Ca (p fj , Pa ,z)-6 Pa% c ia 2 J dxP c {x) z _ £(T _ p2 J{ 2M a )-k*/(2» a y % ~ °' h ^ 

(97) 

where k a = e Q ^(A^ 7 p Q + p^), k' p = e,3 7 (A Q7 p^ + p Q ), a; = p^-p Q , X aP = m a /(m a + mp) = 1 - \p a , a ^ ft. Also 
e Q( 3 = — epa is the antisymmetric symbol with = +1 if (a, 0) is a cyclic ordering of the indices (1, 2, 3). Although 
in Eq. (|97|) formally two Coulomb-modified form factors are present, utmost only one Coulomb-modified form factor 
is needed. In the diagram with i = 0, which describes the neutron transfer (particle 3), the vertex with a = 3 doesn't 
appear. Specifically, the amplitude of the pole diagram describing the neutron transfer 1 + (23) — s- 2 + (13), see Fig. 




FIG. 5: The pole diagram describing the proton transfer in the reaction 2 + (13) — > 3 + (12), which contains the Coulomb- 
modified form factor 1 + 2 — > (12). Notations are the same as in Fig. [3] 
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(23) 




(12) 



+ 



(23) 




(12) 



FIG. 6: The pole diagram describing the nucleus transfer in the reaction 1 + (23) — > 3 + (12), which contains the Coulomb- 
modified form factor 1 + 2 — >• (12). Notations are the same as in Fig. [3] 




FIG. 7: The triangular diagram describing the elastic and inelastic processes in the reaction 1 + (23) — ► 1+ (23). The four-ray 
vertex is the off-shell 1 + 2 Coulomb scattering amplitude. 



[TJ is given by 



R 



,T,-r (1 'i l. , , _ v '<r<r{0)C fm , _^ <r 1 f d f„\ k a L< " k p " X} ( fc a) 

C<3 



z-e°- P l/(2M a )~kl/(2fi a )- 



i = 0, (98) 



where a = 1, (3 = 2 and 7 = 3. Only the vertex /3 + 7 — ► (/S7), which contains nucleus 2, depends on the nucleus 
excitation index a. Correspondingly, the amplitude of the inverse process 2 + (13) -+ 1 + (23), see Fig. [5J is given by 



'-L 



^ii)K7, ( _ p 'trtr(0)£/ / „ . ^_jf_ „<t* 1 [ ka Xq* XcA k a) 



Here, cv = 2, /3 = 1 and 7 = 3. Only the vertex a + 7 — > (07), which contains nucleus 2, depends on the nucleus 
excitation index a. 



In the diagrams with i = 1, 2 (proton or nucleus transfer) the vertex with a = 3 appears only once. Also in 15[ the 
Coulomb interaction in the Coulomb-modified factors and in the four-ray vertex of the triangular diagrams is taken 




FIG. 8: The triangular diagram describing the elastic and inelastic processes in the reaction 2 + (13) — > 2 + (13). The four-ray 
vertex is the off-shell 1 + 2 Coulomb scattering amplitude. 
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(23) 3 (13) 




FIG. 9: The triangular diagram describing the exchange processes in the reaction 1 + (23) — 5- 2 + (13). The four-ray vertex is 
the off-shell 1 + 2 Coulomb scattering amplitude. 



(13) 3 (23) 




FIG. 10: The triangular diagram describing the exchange processes in the reaction 2 + (13) — ► 1 + (23). The four-ray vertex 
is the off-shell 1 + 2 Coulomb scattering amplitude. 

into account in the Born approximation. Here the Coulomb Born amplitude is replaced by the full off-shell Coulomb 
scattering amplitude of particles of the pair 7 = 3. Expression for the Coulomb-modified form factor with the off-shell 
Coulomb scattering amplitude is given above in Eq. (|84]l. 

Now we proceed to the amplitude of the triangular diagram, i — 3. The elastic scattering triangular diagram has 
singularity at forward scattering generated by the off-shell Coulomb scattering amplitude. To eliminate this singularity 
we add and subtract the channel Coulomb scattering potentials in channels a^3, see Eq. (|56l) . The added Coulomb 
potentials can be eliminated by the including Coulomb distorted waves in the initial and final channels a =/= 3. In 
Eq. (A3) [HI the Fourier transform of the screened Born Coulomb potential Vj R \A' a ) should be replaced by the 
unscreened T 7 C , 7 = 3. Then Eq. (A7) 15] for the tr iangular diagrams in Figs [7] and [8] takes the form (a 7^ 3): 

1 

4?£ C2faa (P'<*> P«> z ) = Wo3&r3 {4^:££t>a; z) V 7 %(K) + J dx2 Pc2{x2 ) V c {K) 

-1 

1 

} + J dx2 P ^ ( X2 ) Ftr ift[ tfC 7 ( P " ' Pq 5 Z ) ' 

-1 

t' a < N J '- S " , t a < N J " S " , (100) 
where, owe to /3 = a, p'^ — p' a and t'g =t' a , 

= ^—c°* c a 1 dk k 2 *£«£ g * ( fc )x£f: g w 3 fl011 



T7lJ a S a J a S a GO 

L'J.'^ L a t a f 



/ / \ -1 r> 3 ~=-Ja S a (7(7 

(P ttI V> a ;z) - 16tt j F l ,, , L , 



(p a ;z) 



20 



■DO 



p%/(2M a )-k 2 /(2fx a ) 



\-L' 



Xl y t r '(|A a + k|)|A a + k|- 

X / dxPrAx) — — = o . 

J x z - e° - p'l/{2M a ) - ( A Q + k) 2 /(2 AlQ ) 



(102) 



7-, . J n S a JaS a Cl (7(7 / J \ 

Ftr iZt> a L a t a f (P a ,p a ;z) 



J'„S a a*i 



X J L ±°{k) 



dkk 



x / dxP Cl {x) 





■k|,|A n 



z-*°- p'i/(2M a ) - ( A Q + k)7(2/i Q ) 



*- P 2 /(2Af Q )-fc 2 /(2 Ma ) 



2 T^7 (^71 ^7)^7) 



Here C Q = {S a , L a , J a , t a } and C„ = {S a , L' a , J' a , t' a }, 

A' q =P q -Pq, A q = -A 7/3 A^, x 2 =p a -p«, xi=k-A Q , 



(103) 



(104) 



m p a -m a pp mpp a -m a p p 
^ = — -, k' = Pa+P/3 = Pa+Pfl, k 7 -k; = A^. (105) 

m a /3 m a p 

We use the following notations for the particles in the diagram of Fig. [7] (Fig. [8]): a = a' = 1 (a = a' = 2), 
/3 = P' = 2 (/3 = /J' = 1) and 7 = 3 in both diagrams. The primed particles on the diagrams are the ones after the 
Coulomb scattering described by the four-ray vertex. 

The amplitude (k 7 , k 7 ;z 7 ) is the off-shell Coulomb scattering amplitude of particles j3 and 7 in the triangular 

diagram without the Born term, that is T^{k' , k 7 ; z 7 ) = T9 1 (k' , k 7 ; z 7 ) — VC5 7 (k 7 — k' ), z 7 = z — p T /(2 M 7 ). 
From normalization (|42|) on the energy shell (p a = q a ) we get 



J a S a 



-J a S a aa 

F 

CT=1 *' ,tc=l 



47r E E K:*rll(^;^) = i- ace) 



Then, according to Eq. (f56|) . after applying the two-potential equation the first term 4 n t . £^ ta (p^ ; z) V 7 = C2 (AJ J .) 
in Eq. (|107[> will be replaced by the corresponding Coulomb distorted wave. As a result, we obtain 



1 

p'(3) Ci £2 / cav„/ 



1 

(Pp, p a ;z) = 5 0a 5 a3 5 l3 ^ J dx 2 Pc 2 (x 2 ) {^(A'J F^rf,"^ f° ( p '°>> Pe*;z) - 16 7r 3 F^;^ ta (p Q ;z) 

-1 



+ ^Sttr(P-P- z )}' t'p=t' a <N J '« s «, t a <N J ^. (107) 
Also 

1 

(108) 

The spin-angular part Aclcvf °^ ^ e effective potential is given by Eq. (A2) [Iff, in which q a and q'^ should be 
replaced by p a and p' a . 

The last amplitudes are the exchange triangular diagrams shown in Fig. l9landlT0l 



1 



J„S„ a. 



z- e--pl/{2M a )-kl/{2 i i a ) 



x / d Xl P Cl ( Xl )T-( K ^; Z7 ) ^^^f^ y (109) 

J 1 1 z - e a - p//(2Mr) - (k a + p 2 / 2«« 



z - e° - p//(2 Mp) - (ka + p) 2 /(2 fi a ) 
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We use the following notations for the particles on the diagrams of Fig. [9] (Fig. [T0|) : a = a' = 1 (a = a' = 2), 
P = j3' = 2 (f3 = f3' = 1), 7 = 3, where all primed particles are the ones after Coulomb rescattering. Also 

P' = P/3 + A/3 7 Pa, Xi=k Q -p, 2:2=P^'Pq, 

k 7 = — > A = k 7 - k 7= P" Pa = P/3 - P/3 = PB + P« +P7' ( 110 ) 

P = p/p- The spin-angular part ^■ < clc 2 f^P'/3^P a ^ °^ * ne effective potential is given by Eq. (A8) fl5j |. 

Using the results of Appendix[D]we can simplify the calculations of the effective potentials by combining the neutron 
transfer pole amplitudes with the Coulomb exchange triangular diagram amplitudes taking into account the fact that 
near the pole the exchange triangular diagram has also a pole singularity. Summing up the pole neutron transfer 
amplitude and the pole contribution to the triangular exchange diagram and neglecting the regular at the pole part 
of the triangular exchange diagram, we can replace this two amplitudes by the renormalized neutron transfer pole 
amplitude, that is 

(in) 

where -D^ Q a( *^ is determined in Appendix [D] In Eq. (jllip we took into account that for the pole and triangular 

exchange diagrams Z = Z. The same renormalization procedure can be applied for any tg and t a . Approximation 
(lllip will be used to calculate the angular distributions near the main stripping peak, where the pole neutron transfer 
mechanism gives a dominant contribution, and the results will be compared with the exact approach. 



P — p'j3 — A 7i a Pa, 



K = 



m 8 p' a - m a p' fj 

m n ft 



VI. SUMMARY 



We have derived new generalized Faddeev equations in the AGS form taking into account the target excitations and 
explicitly include the Coulomb interactions. Applying two potential formula we convert the AGS equations to the form, 
in which the matrix elements are sandwiched by the Coulomb distorted waves in the initial and final states. To obtain 
the half-off-shell integral equations we use an off-shell extension. The obtained equations are compact and can be 
solved. We present the final expressions for the modified AGS equations after the angular momentum decomposition. 
We show how to regularize the matrix elements sandwiched by the Coulomb distorted waves. Besides we investigate 
the off-shell Coulomb scattering amplitude in different kinematical regions. We also consider the Coulomb-modified 
form factors and show how to regularize them. After that we investigate the exchange triangular diagram and show 
that its strongest singularity is the pole of the neutron transfer pole diagram. The strongest singularity of the elastic 
scattering triangular diagram with Coulomb four-ray vertex is compensated by the subtracted channel Coulomb 
potential. Thus we have shown that the Coulomb interaction can be taken into account explicitly without Coulomb 
screening procedure. This will allow us to apply the Faddeev formalism for the analysis of the deuteron stripping on 
targets with higher charges, at which the Coulomb screening procedure doesn't work. For NN and nucleon-target 
nuclear interactions we assume the separable potentials what significantly simplifies solution of the AGS equations. 



Appendix A: Partial Coulomb scattering wave function in the momentum space and matrix elements in the 

Coulomb distorted wave representation 



Here we present the expression for the partial Coulomb scattering wave function in the momentum space. This 
wave function has singularity and we demonstrate how this singularity can be regularized when calculating the matrix 
elements, which are given by the sandwiching the diagrams, describing the reaction mechanisms, with the Coulomb 
scattering wave functions in the initial and/or final states. 



First we start from the definition of the Fourier transform of the Coulomb scattering wave function 2l| : 
</£ (+) (p')= lim [ dre- er e- ip '- r ^ {+) (r) 

p v ' e->+0 J 

- lire L(L + iVp) ^ de [(p , _ p)2 + e 2]l-M„ 

= 4 tt J2 V mi (P) Y lmi (pO ^ &/) = ( 2 1 + !) P '(P • P') ^piip')- (Al) 

/ mi I 
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Here rj p is the Coulomb parameter of the interacting particles moving with the relative momentum p, Y\ mi (p) is the 
spherical harmonic function, P/(p • p') is the Legendre polynomial. We can see from Eq. (|A1[) that the Fourier 
transform of the Coulomb scattering wave function is a distribution. 

The expression for the partial Coulomb scattering wave function in the momentum space was found by one of us 
(A.M.M.) [H: 



2 7T In . , iJ,C 



€l(p') = -^e-^T(l + ir, p )e^ Urn 2 i m e -i* "J^"> ^(-1,1 + 1;!-^; 



(p' +p + ie)- 1 +^ P . (p'-p) 2 



(p> -p + ie) 1+ "'p y ' ' /p ' 4pp> 



(A2) 



where 4>f = a? — af , erf is the Coulomb scattering phase shift in the partial wave I, rj p = Z\ Z2 e 2 H12/P is the 
Coulomb parameter for particles 1 and 2 moving with the relative momentum p = £12, 2F\{—1, 1 + 1, ; 1 — 

i rj p ; — ^ ~p} ) is the hypergeometric function, which reduces to a polynomial of order I in the plane z = — ^ . In 

particular, for I = 2 -Fi(0, 1,;1 - ij? p ; - ^p^" ) = 1. 

The function fp p i(p') has singular branching points on the complex plane p' at p 1 = p ± i e and p' = — p ± i e. 
The small imaginary addition ±z e determines the rules for circuiting around the singularities when integrating. If 
the integral containing the Coulomb scattering wave functions is calculated in the analytic form, then no difficulties 
arise, because the presence of the imaginary addend ±i e shifts the singularities from the integration contour to the 
complex plane and lim can be easily taken after carrying out all the integrations. Such a procedure, however, would 

be highly inconvenient in numerical calculations, for in this case it would be necessary to calculate the integrals for 
several continuously decreasing values of e in order to attain a needed accuracy. This procedure, owe to the presence 
of the singularity of the integrand, may become very unstable and even not converging when rj p increases. 

This difficulty can be readily circumvented if, putting e > 0, we regularize the initial integral. Then the result of 
the integration will be stable when integration is performed for lim. The regularization method is taken from [23| . 

e— fO 

To explain the Gel'fand-Shilov method we consider the integral 

1, 

J(A) = lim / dxf(x){x + ie) x , (A3) 

e->+0 J 
a 

where b > and a < 0, ReA = —1, ImA 7^ and /(0) ^ 0. To regularize this integral we assume that —1 < Re A 
and e > 0. Subtracting /(0) from fix) and adding it, we can rewrite integral (|A3|) as 

b 

J{\)=Um o [j dx [f{x) - /(0)] (x + i e )*] + /(0) {b + 1 0)X+1 X ~ { * + t ° )A+1 . (A4) 


J(A) is analytical function of the parameter A in the domain —2 < Re A and A ^ -1. Because J(A) and J(A) coincide 
in the region — 1 < ReA, J(A) is an analytical continuation of J(A) into the domain —2 < ReA. The integrand in 
(|A4[) does not have diverging singularity and the integral can be calculated; the additional term containing /(0) is 
also not singular because A = — 1 + ir\. Eq. (|A4[) can be considered as generalization of the famous equation for the 

b b 

integral containing pole singularity: J dx f{x) /(x — xo — i 0) — P J dx f(x)/(x — xq) + i tt f(xo), where P stands for 

a a 

the Cauchy principal value of the integral. By subtracting and adding df(x)/dx\ x= o x in the integral (|A4[) we can 
continue analytically (|A3|) into the region —3 < ReA. Moreover, in this case the integrand vanishes at x = 0. 
Using the explained regularization we show how it works in practice. As example we consider 

ZI C (PP, P C aQ ) = Jy^ Z ^P^ P) <o(p)- (A5) 


Here subscript denotes / = partial wave. Note that ?p p contains 2Fi(~l,l + 1; 1 — i Vp> ~ ^ V i pp ^ ); which is 
expressed in terms the polynomial of (p — p a ) 2n with < n < I. Hence the terms with n > 1 don't require the 
regularization at the singular point p — p a and the only singular term, which requires regularization, is the one with 
n = 0. That is why it is enough to demonstrate how the regularization works for / = 0. 
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To demonstrate the regularization we consider the effective potential Z(pp, p) given by a simple pole propagator: 

rwt \ 2/zi 2 Zi Z2 e 2 . . 

Z(P "' P) = ~ (P^-P) 2 + ^ ' (A6) 

where p\ 2 is the reduced mass and Z\ Z 2 e 2 is the product of charges of the particles 1 and 2. The partial wave pole 
amplitude at I = is given by 

l 

Z (pp, p) = - p, 12 Z 1 Z 2 e dx- — — -.= <9o(0) A7 ) 

J (P/3-P) 2 + k 2 2p p p 

-1 

£, = {p 2 p + p 2 + K 2 )/(2ppp). From equation [2l[ 

Z sc (pp,p a ) = J ^j%,p)C(P) = -2 m Zi Z 2 e 2 ^ 2 T(l +ir, a ) ^ ~ ■ ( A8 ) 

where 7/ Q = Zi Z 2 e 2 pxilPa is the Coulomb parameter associated with the momentum p Q , which is the on-shell 
momentum. From Eq. (|A8I) using the partial wave expansion |22j | we get 

oo 

Z SC (P0, Pa) = J2 ( 2 Z + !) P '(P^ ' P«) ^(p* Pal)' ( A9 ) 
(=0 



Zr(P0,Pal) 

1 



Pd 

and for / = 



2tt 2 

" a / 2 r(l + ii) a )e < *° Im 



Z,Z + 1; 1 -ir) a ;- 



^PP Pa 



Z { 



SCi 



Pao) 



2_ e -ir Va /2 Y(l + ij7 a )2Im 



^pp - p a + i K 



(A10) 



(All) 



Evidently that analytical expression (IA10[) for Zf (pp, pxj, obtained owe to the very simplified approximation for 
the effective potential Z, can be easily calculated. However, in general case of folding of the partial Coulomb scattering 
wave functions with the effective potentials there are no analytical expressions for Zq [pp, Pao)- I R this case the 
integral must be calculated numerically and regularization of the integrand is required. To show how regularization 
works we compare the results of the numerical calculation of Zq C (pp, p% ) given by the integral representation (|A5[) 
with the analytical expression (| Al 1|) . 

First we present the regularized integral following the Gel'fand-Shilov method described above: 



Zl C {pp lP C a0 ) 
, FR( Pa 



dpp 2 
~2^ 



«+A 



Zoipp, p)^p a0 (p) 



Va Pa 

npp 



Im 



dp 



FR{p) - FR(p a ) - FR'{pa) (p~Pa) 



(p-Pa 



1 



1 



) + ((a + »)'-'«• -(-a + or*-) 



(A + ie) iT >c* (-A + ie)*' 



1 - ir/a 



dpp 2 



2 TV 2 



Zo{pp, p)^ Pa0 (p)- 



(A12) 



Pa+A 



Here we use regularization only in the proximity of the singular point p = p a - To do it we split the integral into 
three terms: the integral from to p a — A, from p a — A to p a + A and from p a + A to infinity. The regularization is 
required only in the integral from p a — A to p a + A, where 



FR(p) = 



Qo(0 

{p+p a y-"i° 



(A13) 
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FIG. 11: (Color online) Comparison of the regularized integrand (dashed red line) in Eq. (|A12[) with the unregularized one 
(solid blue line) as function of the integration momentum p for p < p a , r\ a = 12, k = 0.1 fm - , p a = 1, 1 far - , pp = 1.2 
far 1 and A = 0.25 p a . 
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FIG. 12: (Color online) The same as in Fig. 1111 but for p > p a 



and FR'(p a ) = dFR(p)/dp . Note that here we use regularization procedure subtracting and adding FR(p a ) + 

P=Pa 

FR'(p a ) {p — p a ). In this case the integrand in the regularized integral just vanishes at the singular point. 

To demonstrate how regularization procedure works we performed calculations of analytical Eq. (IA1 1[) and regular- 
ized Eq. (IA12[) for parameters given in the captions to Fig. [TTJ We have chosen the Coulomb parameter r\ a = 12, which 
corresponds to the proton collision with charge 92 (uranium) at the relative kinetic energy sa 1.5 MeV. No Coulomb 
screening procedure would work at such high Coulomb parameter. Results of our calculations are: Zq C (p^, p^q) 
from Eq. (fXTTT) gives -2.009 x 10~ 12 while the regularized Eq. (jM2)l results in -2.005 x 10~ 12 . In Figs [IT] and [T2] 
we demonstrate the behavior of the regularized integrand versus unregularized one. As we can see, the regularization 
completely change the behavior of the integrand making it possible to perform calculations with singular partial wave 
Coulomb scattering wave functions in the momentum space. 

The same procedure can be used to regularize matrix elements with the Coulomb scattering wave functions in the 
initial and final states. Succeeding in regularization of the integrals containing the AGS effective potentials sandwiched 
by the partial Coulomb scattering wave functions, we may conclude that the effective potentials in the AGS equations 
in the Coulomb distorted wave representation can be calculated without using Coulomb screening procedure. It is 
the most important result of our regularization procedure. 
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Appendix B: Off-shell Coulomb scattering amplitude 



In this Appendix using the integral representation of the off-shell Coulomb scattering amplitude we will derive the 
expressions for this amplitude, which can be used in practical calculations in different kinematical regions. We start 
from the standard expression for the off-shell Coulomb scattering amplitude T^(k' a , k a , E a ) [221 l2~il. [25l|: 



T%(k' a ,k a ,E a ) =4nZ p Z 7 



1 



(Bl) 



I{k' a , k a ;E a ) = lim / dxx if >° ^ — ; , (B2) 



>° / ' x(k' Q -k Q ) 2 --^ [^ + 10-^/(2/^)] [£ a + iO-fe2/(2/i Q )] (1- a;) 2 ' 



k Q and k^, are the relative off-shell momenta of particles /? and 7 before and after scattering moving with the relative 

kinetic energy E a , fc" — \jl[L a E a . To distinguish the Coulomb parameter of pair a (between particles (3 and 7) from 
the Coulomb parameter r\ a in the channel a (between particle a and the bound state of the pair (f3"f)) the former is 
denoted by f] a . To underscore that this Coulomb parameter is calculated at fc° we use notation 77° = Zp Z 1 e 2 /_t a /fc°. 
In the case under consideration T^(k' a , k Q , E a ) ^ only for a = 3. 
As we can see T^(k' a , k a , E a ) has the forward singularity at 

A Q = k' a - k a = (B3) 

generated by the Coulomb Born term 4 ir Zp Z 1 e 2 / A 2 . This singularity appearing in the triangular elastic scattering 
amplitude, see Figs [9] and I10| is dangerous when coinciding with the two-body Green function's singularity. Just to 
remove this singularity, we added and subtracted the Coulomb channel potentials in the initial channel a and the final 
channel /?. Applying after that the two-potential equation we obtained the AGS equations in the Coulomb distorted 
wave representation. At A Q — > I{k' a , k a ;E a ) ~ l/|A a |, that is less singular than the Born term. 

It is worth mentioning that the off-shell Coulomb scattering amplitude doesn't have a definite on-shell limit reflecting 
the fact that, owe to the infinite range of the Coulomb interaction, charged particles are not free even when the distance 
between them increases to infinity. Assuming that fc 2 /(2/i Q ) — E a —> we obtain (the on-shell limit in the entry 
channel) 

T C(k> a ,k a ,E a ) ^ r c a {k a , kl)T^ HSH \k' a ,K,E a ), (B4) 
where the so-called Coulomb renormalization factor 

r c a (k a , fc°) = e^/ 2 r(l - ifg) " (B5) 

and the half-off-shell Coulomb scattering amplitude 

T^ HSH Hk' a ,kl,E a ) =4nZ, Z 1 e*e-*^T(l + l f a ) llzi^±l||!L , e . (B6) 

Similar equation takes place when k' a — > fc°, k' a ^ k a . Taking simultaneous limit fc Q , k' a — > fc^ we obtain the 
on-shell Coulomb scattering amplitude. Infinitesimal addend e is required to correctly bypass singularities and for 
regularization of the Coulomb scattering amplitude. In what follows we assume that energy E a has always positive 
imaginary addend i e with e — > 0, that is we replace E a by E a e — E a + ie. 
Let us introduce 



Mi" fc 'a/(2 Ma)] [& a e ~ ^/(2 fl a )] , (B7) 



with b = lim b e , and 



AL = (k^-k Q ) 2 + e 2 . (B8) 
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Then 



T% (k' a , k a , E a ) = lim T Q c e (k' Q , k a , E a , 



where the regularized off-shell Coulomb scattering amplitude is 

1 



A 2 



if] a I t (\s! al k a ;£ ae ) 



(B9) 



(BIO) 



1 

7 e = / dsx**"— , (Bll) 

J x^l^-b^l-x^ 



77° e = Z/3 Z 7 e 2 /i a /y2^ £; Q£ . 

Usage of the regularized Coulomb scattering amplitude allows us to carry out all the calculations with different 
e — > without any problems because the AGS equations are compact at e = (T3, EH- 

Equation (|B10I) is not always the most convenient one for practical applications and below we present alternative 
equations for 7j e , which can be used in practical applications depending on the value of E a . 
(i) E a < 0. 

In this case always b < 0, that is iA^ e -^- (E a - k' 2 J{2 fi a ))(E a - k 2 a /{2 fj, a ))(l - x) 2 > and T° e is regular. Note 
that at E a < 0, owe to the positive imaginary addend, arg(E a +i e) = ir and lim f)® e = Zp Z 1 e 2 fi a /y2 /x Q jE'd e _l7r / 2 



and lim = a; 2 ' 3 ^ e 2 ^ / ^/2^ a \E a \ > . 

(ii) 6 < but £ Q > 0. In this case also no singularities appear in the integrand of the regularized integral over x. 
In general it is more convenient for b < to use an alternative expression for T^ e : 

T ae (k' a , h a ;E ae )=47TZ p Z 7 e 2 -^(E ae -^)(E ae -^) 

1 

„2 



x / dzz*^. IL^J _. ,1512) 

J ^ A 2 , + ^_ (£ ae - k> 2 J(2» a ))(E a€ - fc 2 /(2 Mct )) (1 - x? 



o 



There are no problems in using this explicit equation for T ae (k' a , \z. a ;E ae ) in the Coulomb-modified form factors and 
in the exchange triangular diagram. A strong singularity of lim T ae (k' k Q ; E ae ) in the elastic scattering triangular 
diagram, as has been explained, is compensated by subtracting the channel Coulomb potential. 

(iii) b > (it can be only at E a > 0). In this case the integrand in Eq. (IB11[) has a singularity (zero of the 
denominator): 

x A 2 a -^(E a - ^-){E a - £-) (1 - x) 2 = 0. (B13) 

The roots of this equation are: 

A 2 A 2 / 4b 

^ = 1+ 2? ± 2ff + zf < B14) 

The first root Xi > 1 lies outside of the integration contour over x in (|B11I) at A a > 0, while the second one X2 < 1 lies 
on the integration contour at A Q > 0. The integrand in (|B11|) has a cut connecting the branching point singularities 
at x = and x = oo of the function x lria . Hence we can rewrite /c(k' Q , k Q ; E ae ) as 



k a ;E ae ) = -_- & dxx^>— — — 1 ^— z _ -, (B15) 
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where the integral is taken along the contour C\, which starts at x = 1, encircles the branching singularity at x = and 
ends at x = 1. We assume that the root Xi has a positive infinitesimal imaginary addend, which shifts the root from 
the integration contour to the first quadrant. Now we consider the integral along the closed contour C — C\ +C2 + C3 
shown in Fig. [T51 

Because there are no singularities of the integrand inside and on the integration contour C, according to Caushy's 
theorem the integral taken along C vanishes. Hence, 



1 



1 - e 



-2 7ri)0 



da; ; 



1 



C-2 



■Al 



ae 2 E a 



(Eat 



+ j> dxx lr,a 
c. 



x A 2 



(.E a e-¥t){E ae -^-){i- x y 



(B16) 



The integral over C2 encircles the pole at x = X2 and can be easily taken using Caushy's theorem. The integral over 
contour C3 is taken along the circumference with the radius \x\ = 1. Using substitution x = e 10 we can simplify the 
integral over the contour C3. Then the expression for T Qe (k' Q , k Q ; E ae ) reduces to 



T ae (k' a , k a ;E a< :) = 47rZ (3 Z 7 e 2 



1 



2<, 



A 2 1 - e 



fp 

'la e 



2 7T 



-2 IT f]g, 



A 2 e + 26 £ (l-cos 0) 



4 6 e + Aae) 2 



(B17) 



The integrand in the integral is free of singularities and Eq. (|B17|) can be used to calculate T Q£ (k^, k Q ; E ae ) at b > . 
Note that from this equation we can easily find the on-shell limit for T Qe (k^, k Q ; E ae ) at b —> 



T a t (k' a ,k a ;E ae ) b =° 4 7T Zp Z 1 e 2 -^| Q 

1 — e 



t*5S, Va 2 / 



(B18) 



which coincides with Eq. (IB4|) . 

We have considered different representations of the off-shell Coulomb scattering amplitude and its on-shell limit. 
In the next Appendices we consider how to calculate the coulomb-modified form factors and the triangular diagrams 
using the obtained representations of the off-shell Coulomb scattering amplitudes. 
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Appendix C: Coulomb-modified form factor 



The Coulomb-modified form factor given by Eq. (|84l) consists of two terms. To calculate the integral part we need 
to regularize it. We see immediately that there is a pole in the integrand at 



fc 7 > 0. 



(CI) 



This pole would lead to the logarithmic singularity of the integral when fc 7 = (end-point singularity). However, we 
recall that, according to Appendix [B] fc' 7 — > fc 7 corresponds to the on-shell limit of the Coulomb scattering amplitude 
T^ L (fc 7 , fc 7 ; z 7 ), where it has its own singularity- a branching point. According to Eqs. (|B4[) and (|B5|) . the on-shell 
limit of the partial Coulomb scattering amplitude, see Eq. 



-^7 L~ (^7' ^7> ^7) 



(fc 7 k ^y fl ' 1 T^^ik-y, fc 7 ; z 7 ), 



(C2) 



where ?) 7 = Z Q e 2 /i 7 / ^/2 /i 7 z 7 . At z 7 



i? 7e = E^ + ie, where e — > 0, ?) 7 = fy" . Note that in practical calculations 



£ + -p 2 /(2M 7 ), 



£* + = E + iO. Then the integrand in the integral part of Eq. (f84|) has a singular behavior 



■^yi,^' ■ 2 "y) fe' 
z 7 -F 2 /(2M 7 ) 



7 , fc 7 , z 7 ) 



z 7 -/c' 2 /(2 M7 )]!-»h 



(C3) 



Thus the integrand in the Coulomb-modified form factor (|84p has a branching point singularity rather then a simple 
pole. As it has been demonstrated in Appendix this type of singularity can be easily regularized using Gel'fand- 
Shilov method [23| . We apply here this technique to regularize the integral in the Coulomb- modified form factor. 

First off all such regularization is required only at Re£ 7 > because at Rez 7 < and Imz 7 =0 z 7 — fc'^ < 0, 
that is the integrand is regular function. To regularize the integral at Re z 1 > and Im z 7 — > +0 we rewrite the 



Coulomb-modified form factor 



in the form 





dk' k' 2 ^ L ^ kl1 ^' Z ^ *i ^7) 



,/ 2 ^^(^j^^tOXl^ (^7) 
7 7 [z 7 -fc' 2 /(2 M7 )]i-^ 



7 7 [z 7 -fc' 2 /(2 Al7 )]i-^ 



(CM) 



In the first integral arg(z 7 — fc 7 7 /(2 /z 7 )) — > if Im z. 
Then we can rewrite Eq. (|C4I) as 



dtftl "fa) = xtft/fa) + (5 73 ^ 



0. In the second integral z 7 —fc' 7 /(2 /i 7 ) = e i7r [fc' 7 /(2 /z 7 ) — z 7 ] 



7 7 [z 7 -fc' 2 /(2 Ai7 )]i-^ 



„, ,,2 ^lA^"i->^>^i)Xl~Z (^7) 

dfc., fc^, 1 — 1 e 

7 7 [fc' 7 /(2 M7 )-z 7 ]i-^ 



2fe, 



j, / , / 2 ^L T (^7' ^7' Xl^, (^7) 

QrC.. rb., o " 

7 7 [fc' 7 /(2 M7 )-z 7 ]i-»H 



(C5) 



Assuming that Re i fj-y > we can rewrite Eq. (|C5|) in the form, in which the singularity is weakened enough to 
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make it integrable at Rei 7/ 7 — > 0: 



r dk' k! ^ "^T^-r^T" ^7' ^) *7 ^TlJ^'^^l) Xi^t^ (^7) 



7 7 



[z 7 -fc' 7 /(2 Al7 )]i-% 



+ f^-z^ 1 " 1 k 1 T^ L ^{k 1 ,k 1 -,z 1 )xLt vh) 



, dfc' fc' ^ ^t^t^ 7 ' ^7)Xi 7 t 7 (& 7 ) ^ 7 T^^^, fc 7 ; £ 7 ) (^7) 



7 7 



[fc' 7 /(2/x 7 )-i 7 ]i-i^ 



+ i^e-^ (3z 7 r^ fc 7 f 7 c L7 (fc 7 ,fc 7 ;z 7 ) X ^t/(fc7) - e" 7 ^ / dfc 7 fe 



iry 7 



./ / 2 ^L 7 (^7' ^7! ^7) X.L 7 ^ (^ 7 ) 



[fc' 7 /(2 M7 )-z 7 p-^ 



(C6) 



Eq. (|C6[) can be used to calculate the Coulomb- modified form factor. We have shown how to regularize the Coulomb- 
modified form factor in the proximity of the singularity fc 7 = fc 7 . This regularization allows one to calculate the 
Coulomb-modified form factor. 

Now we consider a special case of a simple separable form factor 



A—.. 



(C7) 



for which we will derive analytical expression for the on-shell limit of the Coulomb-modified form factor. The Coulomb- 
modified form factor in this case can be written in the integral form [26| : 



J-y5-. (J I 



-4fc 2 



dx ■ 



\x B + - a - a 



TOT- (C8) 



Here, 



B 



tilt' -iky' 



(C9) 



From Eq. (|C8I) clear that the Coulomb-modified form factor has singularity in the on-shell limit fc 7 — > fc 7 corre- 
sponding to a -> 0. For the separable form factor (|C7[) at L 7 = the on-shell limit of the Coulomb-modified form 
factor was found in [l8| (see Appendix C): 



9ot y 



2 fj-y arctan fc 7 //3 Q 7 7 CT / ^'7 ^'7 

l_ e -2^ 7 e T V 2 fc 7 



(CIO) 



Here, we extend this expression for arbitrary L 7 . To do it we note that Eq. (|C8[) can be rewritten as 



-4fc 2 



L 7 + l 1 



dx ■ 



«2 ~\~ L ^ 



where 2^2 are the roots of equation x B + — a — a -1 = what is equivalent to a; 2 — -g (a + a ) x + -gy = 0: 

(C12) 



#1,2 — — — ± 



2B 



2B 



Evidently that at k that is a — > 0, 



a->0 1 a->0 a „ 

xi = — - } oo , x 2 = — -> 0. 
ai3 i( 



(C13) 
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Also x lf, ~<~ 1 in the integrand of Eq. (|C11[) has a branching point singularity at x=0 and oo, which are connected 
by a cut. Then we can apply the method described in Appendix [B] First we transform the integral in (|CJ11|) to the 
integral over contour C±, in which the integration contour is taken along the contour starting at x = 1, encircling the 
branching singularity at x = and ending at x = 1, see Fig. Q21 After that we consider the integral along the closed 
contour C = C\ + Ci + C3. Because there are no singularities of the integrand inside and on the integration contour 
C, according to Caushy's theorem, the integral taken along C vanishes. Hence 



dx 



x 



1 



[(x — Xl)(x — X2)] Ll+1 



1 — g-2 7rr) 7 



dx 



[(X ~ X\)(x — X2)] Ll 



+ d> dx 



c 3 



[(x — Xl)(x — X2)] L ~* 



(C14) 



We assume that the root Xi has a positive infinitesimal imaginary addend, which shifts the root from the integration 
contour to the first quadrant. In the integral over contour C2 we rewrite 



lim 



1 



lim 



1 



[(x — xi)(x — X2 — e)] L -t +1 L 7 ! (x — xi) L ~i +1 e-todf i i x — X2 — e 



(C15) 



We take into account that at a — > X2 — > while x\ — > 00, that is at a — » contour C2 encircles the pole at x ~ X2 + e 
and can be reduced to the residue in the pole: 



1 



da; ■ 



x 



.1 ?}-y ~\~ Ly 



1 



1 — e 27rf >i J ' [(x — Xi)(x — x 2 )] L ^ +1 
c 2 



1 ,. d L '< 
hm 



1 - e -2"h LI e ^ de L 



dx ■ 



r i fjy+L-, 



Co 



[{x — xi)] L i +1 (x — X2 — e) 



2tti 



1 d L ^ (x 2 + e) l ^+ L -< a ^0 

hm — 



2?r i 



1 - e - 2m "^ L 7 ! £^0 de L ~> (x 2 ~x 1 + e) L T+i l - e -^fh L 7 ! 



I 1 

( 1) L ^ + 1 a if l~>+ L -> + 1 g-ifii+L^ + l TT / 



irjy + n), 
(C16) 



L T =0 

where ]~[ (i?) 7 + n) = % fL. Also at a — > a = (fc 7 — k 7 )/(2 fc 7 ) and i? = (/3 L 7 4 7 + i fc 7 )/(/3 L 7 f ** — i fc 7 ). 

n=0 



Taking into account all the factors in front of the integral in Eq. (jClip 

1 1 1 



i 7}^ 1 



-4fc 2 



dx ■ 



x 



2tt 



B L ~t +1 1 - e" 2 ^ J [(x - xi)(x - x 2 )] L ~< +1 
c 2 

L-, 



1 — e~ 2vf i-t 



V~f 2 ?7-y arctan / /3 ' L 



v 2L / 



K\ if i-i 1 1 



L 7 ! i?7a 



n (^7 + ™)xl^' t (>7)- 



(C17) 



71=0 



Let us consider the integral over C3. In this integral the contour goes along the circumference with |x| = 1. Hence, 
along C3 x — e %v and dx = i d <p e l v . 



1 



dx ■ 



x 



2 71 



a->0 



— (-aS)^ 1 / d^e 



1 — e 2nf i~t / [(x — xi)(x — X2)] L ' I+1 1 — e 

C 3 

= i(-l) L i (aB) L -< +1 —. 

Taking into account the factors in front of the integral in Eq. (jCllj) we get 

1 



-V-t v 



i fjry k. 



-4fc 2 



dx ■ 



i fj 1 +L 1 



l_ e -2Kr ) -,J [(x - Xi)(x - X2)]^ +1 

C 3 



(ci8) 



(C19) 
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Then we arrive at the final equation of this Appendix defining the on-shell limit of the Coulomb-modified form 
factor for the separable form factor (|CJ7|) : 



kry—tk-y Iht)^ 2-fjry arctan fe-y//^ 7 

I _ g-2 TTfly 



-±k 2 



L-, + 1 I 



B L ->+ 1 



dx 



x 



[(x — Xi)(x — X2)] L ~< + 1 



V 2fc 7 ) 



' kl M ^ ' ' n^+^X^t/CM, k 7 >0. (C20) 



LJ. ifjx 



n=0 



k^ . 



This 



Thus we have shown that the Coulomb-modified form factor has a branching point singularity at fc 7 
Coulomb-modified form factor is needed to calculate the amplitudes of the diagrams describing the proton and 
nucleus transfer, see Eq. (j9Tj) . Assume that the Coulomb-modifed form factor is g° (k a ), which has the on-shell 

singularity at k a — k°, where k° = y7/^ (z a — eF), z a = z — p 2 a / (2M a ). Then we can rewrite Eq. (|97|) as 



R 



Cl (Pff'Pcxi z ) = R Q,'(« >L '(p'i3'P°'' Z ) = S Pa c c/ c C= 9 / dX 



■(i)C 



K La k^xi;{k' p )~gl{k a 
[za-eP-ki/^)} 1 -^ 



* = 0, 1, 2, 

(C21) 



where g° (k a ) is regular at k a — > k a a . Thus the integrand, owe to the presence of the Coulomb- modified form factor, 
has a branchin g p oint singularity rather than a pole and the integral can be easily regularized using the Gel'fand- 
Shilov method [23| applied earlier. Note that, although derivation of Eq. (|C20j) has relied on the explicit form (|C7p 

of the separable form factor \l t^ik-y), the result is nevertheless valid for arbitrary nonsingular at the origin form 
factor since any such form factor can be represented as linear combination of functions of type (|C7|) . Thus we have 
demonstrated that the Coulomb-modified form factor does not create any problems in practical calculations. 

Finally, it is worth mentioning that the on-shell limit of the Coulomb-modified form factor at k y = for L 7 = was 



found in [13: g^" (fc,) *' 



Similar consideration can be done for L 



k 2 , that is the Coulomb-modified form factor vanishes at the on-shell limit at fc 7 = 0. 



> 0. 



Appendix D: Pole singularity of the triangular exchange diagram and Coulomb renormalization of its 

strength 

In the previous section we considered the Coulomb-modified form factor and how the presence of the off-shell 
Coulomb scattering amplitude affects it. The off-shell Coulomb scattering amplitude also is needed to calculate the 
direct and exchange triangular diagrams. The discussion of the direct triangular diagram has been done in Section 
Ivl Here we consider the exchange triangular diagrams shown in Figs [9] and 1101 

In this section we show that there is a pole singularity of the exchange triangular diagram, which allows us to 
rewrite the amplitude of this diagram as the renormalized pole neutron transfer diagram plus the nonsingular at the 
pole term. To show it let us consider, for example, the amplitude of the exchange triangular diagram shown in Fig. 
IH1 For simplicity, we neglect the orbital momenta in the three-ray vertices and the spins, although the final result, 
which we present below, is valid for the general case. The amplitude of the exchange triangular diagram is given by 

where a = 1, (3 — 2 and 7 = 3; k 2 a — 2 [i a z a and z a ~ z — p^/ (2 M a ); = 2 [ip zp and zp = z — p' 2 p/(2 Mp); 
p 7 + p a and p 7 + -^ J - p'p are the relative momenta of particles in the three-ray vertices (^7) — > (3 + 7 and 

(aj) — > a + 7 in the triangular diagram. Also k 7 = p Q + (m a /m a p) p 7 (k 7 = ~p'p — {mp/m a p) p 7 ) is the relative 
momentum of particles a and /3 on the diagram before (after) the Coulomb scattering with the transfer momentum 
in the four-ray vertex = k 7 — k 7 = p a + p^ 4- p 7 (see also Eq. (|110[) V 

The closest to the physical region and the strongest singularity of the triangular diagram is the one generated by 
the coincidence of the singularities of the propagators — (p 7 + p a ) 2 — and fc^ — (p 7 + ■^ J - p'p) 2 = and the 
forward singularity of the off-shell Coulomb singularity A^ = of the Coulomb scattering amplitude. To show how 
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these singularity of the exchange triangular diagram appears and to simplify consideration we replace the off-shell 
Coulomb scattering amplitude T^(k^, k 7 ; z 7 ) of particles a and /3 by the Born Coulomb amplitude, which is the 

Fourier transform of the Coulomb potential 4 7T/A' 2 . Then the amplitude of triangular exchange diagram simplifies 
to 

' 1 / (2ir) 3 h 2 - fr> A- -22- n' "l 2 A' k 2 — (rt I r» \ 2 

J (2 ») 3 k% - (AJ, - kj)" A* fcj - (AJj - k,,) 2 

where we used the substitution p 7 = — p Q — p«. Also k Q = + (m^/m^ 7 ) p a is the relative momentum of 
particles /3 and 7 in the three-ray vertex (^7) — > /3 + 7 in the diagram of Fig. [1] with a = 1, /3 = 2 and 7 = 3. 
Similarly = p a + {m a / m ai ) p^ is the relative momentum of particles a and 7 in the three-ray vertex (cry) — > a + 7 
of the same diagram. Now we rewrite 

k 2 a - (A' a - k Q ) 2 = a a + 2 • k Q - A'l = a a [l + 2 1 • k Q - a a t 2 ] , (D3) 

where we introduced a a = k a — k 2 and used the substitution 

K = <r a t. (D4) 

Similarly 

kp — (A„ — k^,) 2 = a a [— +2t-k' - a a t 2 ]. (D5) 

f^a 

Here we took into account that from the energy-momentum conservation in both three-ray vertices of the diagram in 
Fig. Q] follows 

<Tp = —(T a . (D6) 
Ma 

Because we consider the singularity of the exchange triangular diagram generated by the coincidence of zeroes of 
three denominators (pinch-point singularity) in Eq. (ID2[) . we use the substitution (|D4[) obtaining in the leading order 

Z^( PpiPa ) = -4nZ ? Z a e [ \J ^-^-^ (D7 ) 

Thus we have shown that the strongest singularity of the amplitude of the exchange triangular diagram is a pole 
singularity at 

a a = 0. (D8) 
The same singularity has the pole diagram in Fig. [T] Moreover 

z^( P ' , Pa ) ~ r *S(- k g)x°(- k °) i (D9) 

L a a i 

in Eq. (|D7j) is the amplitude of the pole diagram (we neglect the spins and angular momenta). We can conclude 
from the simple consideration presented here that near the singularity (|D8I) the amplitude of the exchange triangular 
diagram behaves as renormalized amplitude of the pole diagram in Fig. [TJ 

Z^(p' ,p a ) ° a =°DZ(°\p' p ,p a ) + z£l(pp,p a ), (D10) 

where D is the renormalization factor determining the strength of the pole singularity. The additional term 
Zrel(p'p,p a ) is regular at a a — 0. 

A general expression for the renormalization factor for the exchange triangular diagram containing the full T 7 
Coulomb scattering amplitude rather than the Born Coulomb amplitude was obtained in (27l - [29l |. Summing up the 
neutron transfer pole diagram and the corresponding exchange triangular diagram we obtain the renormalized pole 
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diagram plus the additional term from the exchange triangular diagram, which is regular at a a = 0. For the diagram 
in Fig. [5] the renormalization factor, which determines the strength of the pole singularity of the triangular exchange 
diagram, is 



T) a — — 1 



y/-m^m a (z a - e CT ) + y 1 -m ai mp(zfi - e a ) - iJ m 7 m0 a E^ (/jaC 



(Dll) 



rp(T rp<J 



M 



2 mf 



2 m r 



(D12) 



(D13) 



(9C 3 "Co 



(D14) 



2 A/„ 



z = z — , and M. = m a + mp + m 7 . For the diagram in Fig. 1101 



IT" 7 

"Co 



(D15) 



Thus we have shown that the off-shell Coulomb scattering amplitude doesn't cause any problem in calculation of 
the exchange triangular diagram. 
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